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In this paper we study the scattering theory associated with the pseudofermion dynamical theory 
for the Hubbard chain. While for the electronic basis the problem is non perturbative and strongly 
correlated, in terms of pseudofermions the spectral properties are controlled by zero-momentum 
forward scattering only. Indeed, we find that each ground-state - excited-energy-eigenstate tran- 
sition corresponds to a well defined set of elementary two-pseudofermion zero-momentum forward- 
scattering events. An important point of the theory is that independent 77-spin 1/2 holons and 
1} ' spin 1/2 spinons are neither scatterers nor scattering centers. Instead, the scatterers and scat- 

5^ , tering centers are spin-less and 77-spin-less cO pseudofermions, 77-spin-zero 2^-holon composite cv 

^ ■ pseudofermions, spin-zero 2^-spinon composite sv pseudofermions such that v — 1,2,3, and the 

corresponding pseudofermion holes. Similarly to chromodynamics, where all quark-composite physi- 
cal particles are color-neutral, for the pseudofermion dynamical theory all 2^-holon (and 2^-spinon) 
composite pseudofermion scatterers and scattering centers are 7/-spin- neutral (and spin-neutral). 
Thus, the pseudofermion S matrix is a mere phase factor, which is behind the simple form of the 
\ pseudofermion anti-commutators and the simplification of the study of the finite-energy spectral 

and dynamical properties. The pseudofermion S matrix is expressed as a commutative product of S 
matrices, each corresponding to an elementary two-pseudofermion scattering event. This commuta- 
tive factorization is stronger than the usual factorization associated with Yang-Baxter Equation for 
the original spin 1/2 electron bare S matrix. Our results reveal the scattering mechanisms which 
control the exotic finite-energy spectral properties of the low-dimensional complex materials and 
correlated systems of cold fermionic atoms on an optical lattice. Importantly, the exotic scatter- 
ers and scattering centers predicted by the theory were observed by angle-resolved photoelectron 
spectroscopy in low-dimensional organic metals. 

^ . PACS numbers: 70 

O 

O '■ !• INTRODUCTION 



The one-dimensional (ID) Hubbard model is solvable by coordinate Bethe ansatz (BA) |l|, |2J. The solution of 
the model can also be achieved by the inverse-scattering algebraic BA Jj|. In terms of electronic scattering the 
model describes a very complex non-perturbative strongly correlated problem. The main goal of this paper is to 
show that in terms of the pseudofermions associated with the holons and spinons introduced in Refs. Q |^ @ the 
scattering problem considerably simplifies and involves two-pseudofermion zero-momentum forward scattering only. 
O ■ In the last twenty years the low-energy behavior of the model correlation functions has been the subject of many 
> I studies iTtlaJplHi [ill [ij, [IJ . Recently, the use of the pseudofermion representation extended such studies to finite 
j^j ■ energy EF 

The generalization of the electron - rotated-electron unitary transformation introduced in Ref. [13 for all values of 
the on-site repulsion U plays an important role in the construction of the holon and spinon representation of Ref. Q . 
Such a representation is faithful for the whole Hilbert space, including the subspace generated by application onto 
the BA-solution states of the off-diagonal generators of the spin and 77-spin SU(2) algebras In turn, there 

is one pseudofermion representation for each initial ground state whose subspace is spanned by the excited energy 
eigenstates contained in the excitations generated by application of one- and two-electron operators onto the former 
state 0, ^3 . The pseudofermion dynamical theory (PDT) 0, 0, is a suitable starting point for the study of 
the dynamical and spectral properties of the model for all values of the momentum, energy, on-site repulsion U, and 
electronic density n. However, the relation of the PDT to the scattering mechanisms remains an open question. The 
theory is a generalization for all values of the on-site repulsion U of the U/t >> 1 method of Ref. j2jj. Here t is the 
first-neighbor transfer integral. 

The finite-energy PDT reproduces the well-known behavior of spectral and correlation functions in the limit of 
low energy , which was previously obtained 0, H, 0, El El 13 E3 by use of methods such as conformal-field 
theory [2IJ and bosonization • The theory was successfully applied to the description of the unusual finite-energy 
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spectral properties of low-dimensional complex materials [23i, |2J, |25j : For the one-electron removal spectral function 
the singular spectral features predicted by the PDT show quantitative agreement for the whole energy band width 
with the peak dispersions observed by angle-resolved photoelectron spectroscopy in the quasi- ID organic conductor 
TTF-TCNQ (Results for the TTF-TCNQ spectrum consistent with those of the PDT were obtained by the 

dynamical density matrix renormalization group method |26j.1 Moreover, the theory was also used in the description 
of the phase diagram of other low-dimensional complex materials [25| and is of interest for the study of the spectral 
properties of the new quantum systems described by cold fermionic atoms on an optical lattice |27l | . (New experiments 
involving cold fermionic atoms [such as 6 Li] on an optical lattice formed by interfering laser fields are in progress 

We are able to calculate explicitly the pseudofermion and pseudofermion hole S matrices and overall phase shifts 
and to clarify how these quantities control the finite-energy spectral properties. Interestingly, the unusual independent 
charge and spin spectral features observed by angle-resolved photoelectron spectroscopy in low-dimensional organic 
metals |23l \24\ correspond to the exotic independent charge and spin pseudofermion scatters and scattering centers 
introduced here. This paper contains a detailed presentation of the preliminary results presented in short form 
elsewhere (29j. 

The relation of the holon, spinon, and pseudofermion description of Refs. II , 1 5 IS IfH ITU ITi| used here to the 
conventional holon and spinon representation and scattering theory of Refs. 30, 31, 32] was very recently clarified in 
Ref. |33j | . Such an investigation confirms that both representations are faithful and thus that there is no inconsistency 
between the two corresponding definitions ofquantum objects. Moreover, that study confirms that the holon, spinon, 
and pseudofermion description of Refs. 0,00, Q, 0, E| is the most suitable for the study of the finite-energy 
spectral and dynamical properties. 

The paper is organized as follows: In Sec. II we introduce the ID Hubbard model and summarize the basic 
information about the pseudofermion description needed for our studies. The general pseudofermion scattering theory 
and the pseudofermion S matrix are introduced in Sec. III. In that section we also discuss the relation between the 
pseudofermion scattering theory and the spectral and dynamical properties. In Sec. IV we introduce and study the 
pseudofermion phase shifts. Finally, Sec. V contains the concluding remarks. 

II. THE ID HUBBARD MODEL AND THE PSEUDOFERMION DESCRIPTION 

In this section we introduce the ID Hubbard model and summarize the concepts and results concerning rotated 
electrons and the pseudofermion description 0] that are needed for our studies. 

A. THE ID HUBBARD MODEL AND ROTATED ELECTRONS 

The exotic quantum objects associated with the BA solution of the ID Hubbard model are related to the electrons 
through the rotated electrons 0, ■ Let us start by introducing the model. 

1. THE ID HUBBARD MODEL 
In a chemical potential fj, and magnetic field H the ID Hubbard Hamiltonian can be written as, 

H = H S 0(4) + ]T Ma S z a , (1) 



where 



Hso(4) = Hh 



2 L 2 



H H = f + UD. (2) 



Here Hh is the "simple" Hubbard model, T = —t Ylo-=l 1 Sj=i c j a c i+ 1 . a + ^- c - * s ^ ne kinetic- energy operator, D = 
12f=i c ] t c i> T c ] x C J' J- = 5Zj=i T 1 i s the electron double-occupation operator, and the operator hj t a = a Cj, a 
counts the number of spin-projection a electrons at lattice site j. The operator a (and Cj, a ) that appears in the 
above equations creates (and annihilates) a spin-projection a electron at lattice site j = 1, 2, ...,N a - We consider that 
the number of lattice sites N a is large and even. On the right-hand side of Eq. QJ, /j, c = 2/x, ^ s = 2[ioH, is the 
Bohr magneton, and the number operators S* — — \[N a — N] and S% = — |[iVj — N{\ are the diagonal generators of 
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the 77-spin and spin SU(2) algebras 0,0, respectively. Here the electronic number operators read N = J2a=*\ 1 
and N a — ™j,<r- The momentum operator is given by P — X) CT =t 1 Sfe ?v(&) k, where the spin-projection a 

momentum distribution operator reads n a (k) — c\ a c^ a and the operator c\ a (and Ck,a) creates (and annihilates) 
a spin-projection a electron of momentum k. 

Throughout this paper we use units of both Planck constant h and lattice constant a one. We denote the electronic 
charge by — e, the lattice length by L — N a a = N a , and the 77-spin value r\ (and spin value S) and 77-spin projection 
rj z (and spin projection S z ) of the energy eigenstates by S c and (and S s and Sf), respectively. The Hamiltonian 
Hso(4) given in Eq. (0 commutes with the six generators of the 77-spin and spin SU(2) algebras and has £0(4) 
symmetry |0, ^| . While the expressions of the two corresponding diagonal generators were given above, the off- 
diagonal generators of these two SU(2) algebras read §1 = l}* c\ , c\ T and S c — c j, T c i, I f° r ^ s P m 
and 5j = J^j °j 1 c j, T anc ^ ^ = c ] i c j, I f° r spin. The BA solvability of the ID Hubbard model is restricted 
to the Hilbert subspace spanned by the lowest-weig ht states (LWSs) H or highest-weight states (HWSs) [3j of 
the 77-spin and spin algebras, that is by the states whose S a and S* numbers are such that S a = —S* or S a = S*, 
respectively, where a = c for 77-spin and a = s for spin. Such states have electronic densities n — N/L and spin 
densities m — [N-\ — N\\/L in the domains < n < 1 and < m < n, respectively. The description of the states 
corresponding to the extended domains 0<ti<1 ; 1 < n < 2 and — n < m < n ; — (2 — n) < m < (2 — n), 
respectively, is achieved by application onto the latter states of off-diagonal generators of the 77-spin and spin SU (2) 
algebras [J. The scattering processes studied in this paper result from ground-state - excited-state transitions. For 
simplicity, here we consider initial ground states with densities in the domains < n < 1 and < m < n, respectively. 
(Some of our results correspond to initial ground states with densities in the ranges < n < 1 and < m < n.) 

2. ROTATED ELECTRONS 

Each lattice site j = 1, 2, N a of the model (JTJ can either be doubly occupied, unoccupied, or singly occupied by a 
spin-down or spin-up electron. The maximum number of electrons is 2N a and corresponds to density n — 2. Besides 
the TV electrons, it is useful to consider N h = [2N a — N] electronic holes. (Here we use the designation electronic hole 
instead of hole, in order to distinguish this type of hole from the pseudofermion hole.) Our definition of electronic 
hole is such that when a lattice site is unoccupied, we say that it is occupied by two electronic holes. If a lattice site 
is singly occupied, we say that it is occupied by an electron and an electronic hole. If a lattice site is doubly occupied, 
it is unoccupied by electronic holes. The same definition applies to the rotated-electronic holes. 

The electron - rotated-electron unitary transformation maps the electrons onto rotated electrons such that rotated- 
electron double occupation, non-occupation, and spin-up and spin-down single occupation are good quantum numbers 
for all values of U/t 0. The lattice occupied by rotated electrons is identical to the original electronic lattice. 
We call Cj a the electrons that occur in the ID Hubbard model (|TJ and J2J), while the operator ct a such that 
Cj a = V^(U/t)Cj a V(U/t) represents the rotated electrons, where V(U/t) denotes the electron - rotated-electron 
unitary operator. Similarly, Cj a = V(U/t) c] a V^(U/t). Note that for m = cj g and cj a are only identical in the 
U/t — > 00 limit where electron double occupation becomes a good quantum number. The operators V'(U/t) and 
V(U/t) are uniquely defined for all values of U/t by Eqs. (21)-(23) of Ref. Q. The electron - rotated-electron unitary 
transformation was introduced in Ref. 0|. The rotated-electron double occupation operator D given in Eq. (20) of 
Ref. Q commutes with the ID Hubbard model. We denote the rotated-electron double occupation by D r . 

B. THE PSEUDOFERMION DESCRIPTION 

Here we summarize the pseudofermion properties that are needed for the studies of this paper. The pseudoparticles 
studied in Refs. 0,0 and the pseudofermions used in the investigations of Refs. [THllTfil | are closely related. While the 
pseudoparticles have discrete bare-momentum values qj such that q-j+i —qj = 2ir/L, the corresponding pseudofermions 
have canonical-momentum values q~j = qj +Q% v {<lj)/L. (The designation "bare-momentum" follows from the discrete 
values q 3 in units of 2ir/L corresponding to quantum numbers, as given in Eq. (B.l) of Ref. 0; The designation 
" canonical- momentum" stems from the analogy of the momentum shift Q^uilj) I 'L with the shift within the canonical 
momentum of electrons in the presence of a vector potential.) Here Q% v {qj) / L is the functional given in Eq. (14) of 
Ref. and av labels the pseudofermion branch, as discussed below. Although that functional is of the order l/L, 
the discrete canonical-momentum are such that (jj+i — q~j — 2tt/L + 0(1/L 2 ). Except for the slightly different discrete 
canonical-momentum values q~j and discrete bare- momentum values qj, the pseudofermions have the same properties 
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as the corresponding pseudoparticles. For instance, they have the same values of charge, 77-spin, or spin and for the 
branches other than the cO branch, also the same holon or spinon contents. The pscudofcrmion description refers to a 
Hilbert subspace called in Ref. @ pseudoferrmon subspace (PS). (All one-,two-, and any other finitc-number-electron 
excitations are contained in the PS.) In the PS the energy eigenstates are described by the same pseudoparticle Q 
and pseudofermion [f| occupancy configurations and the av pseudoparticles and av pseudofermions are related by a 
unitary transformation 0. Thus, the basic pseudofermion properties summarized below have many similarities with 
the corresponding pseudoparticle properties studied in Refs. 



1. cO PSEUDOFERMIONS, COMPOSITE PSEUDOFERMIONS, YANG HOLONS, AND HL SPINONS 

A key result needed for our study is that the energy eigenstates that span the PS can be described in terms of 
occupancy configurations of holons, spinons, and cO pseudofermions jj,|6|. We recall that the holons and spinons 
considered here are different from those of the conventional spinon-holon representation used in the studies of Refs. 
|30l 1-3 lL l32| and that the relation between the two alternative holon and spinon representations is clarified in Ref. 
[33|. For the simplest excited energy eigenstates, the holon (and spinon) of the conventional representation involves 
mixing of the cO pseudofermion hole and Yang holon (and si pseudofermion hole and HL spinon) considered below. 
The holons (and spinons) introduced in Ref. [4( have 77-spin 1/2, 77-spin projection ±1/2, charge ±2e, and spin zero 
(spin 1/2, spin projection ±1/2, and no charge degrees of freedom). We use the notation ±1/2 holons (and ±1/2 
spinons) according to the value of 77-spin projection (and spin projection). The rotated-electron double occupation D r 
equals the number of —1/2 holons. Within the description of charge transport in terms of electrons (and electronic 
holes), the cO pseudofermions carry charge — e (and +e) and have no spin or 77-spin degrees of freedom. Moreover, 
the cv pseudofermions (and sv pseudofermions) are 77-spin zero (and spin zero) composite objects of an equal number 
v = 1, 2, ... of —1/2 holons and +1/2 holons (and —1/2 spinons and +1/2 spinons). Within the description of charge 
transport in terms of electrons (and electronic holes), the cv pseudofermions carry charge —2ve (and +2ve) where 
v = 1, 2, .... In this paper we use the notation av pseudofermion, where a = c, s and v = 0, 1, 2, ... for the cv branches 
and v = 1,2, ... for the sv branches. The ±1/2 holons (and ±1/2 spinons) which are not part of 2z/-holon composite cv 
pseudofermions (and 2^-spinon composite sv pseudofermions) are called ±1/2 Yang holons (and ±1/2 HL spinons). 
In the designations HL spinon and Yan g ho lon, HL stands for Heilmann and Lieb and Yang refers to C. N. Yang, 
respectively, who are the authors of Refs. |l8Lll3 |. We denote the number of av pseudofermions by N av and the number 
±1/2 Yang holons (a = c) and ±1/2 HL spinons (a — s) by L ai ±i/ 2 . Note that Y c o equals the number of rotated- 
electron singly occupied sites, [iV a — iV c o] equals the number of rotated-electron doubly occupied plus unoccupied sites, 
and L Qj ±i/2 = S a + S^. We call A/ Qi ±i/ 2 the number of ±1/2 holons (a = c) and ±1/2 spinons (a — s) such that 
M Q , ±1/2 = L at±1/2 + Y^ivNau. These numbers are given by M c ,_ 1/2 = [N — N c0 }/2, M c , +1/2 = [N h - Y c0 ]/2, 
M s ,_ 1/2 - [A^o-A^+Yj/2, and M. t+x /2 = [Y c0 + Y T - Y x ]/2. Furthermore, M a = [M a< _i/ 2 + M Q , +1/2 ] denotes the 
number of holons (a = c) or spinons (a = s) such that M c = [N a — N c0 ] and M s = N c0 and L a = \L a< _xli + L a , +1/2] 
denotes the number of Yang holons (a — c) or HL spinons (a — s) such that L c = 2S C = 2r\ and L s — 2S S — 25. 
An important point is that for the ground state and densities such that < n < 1 and < m < n one finds that 
N c0 = N, N sl = N u M Cj+1/2 = L C:+1/2 = [N a - N], M Sj _ 1/2 = N u M Si+1/2 = Y T , L s , +1/2 = [jVj - N^, and 
N av = M Ct _x/2 = L a _i/ 2 — for av ^ cO, si and a = c, s. 

Often in this paper we use the notation av ^ cO, si branches, which refers to all av branches except the cO and 
si branches. Moreover, the summations (and products) E al /=cO, s i> and E a ^ c o, s i ( and U. av , Uau=cO,si> and 

riai/^co si) run over ai l av branches with finite av pseudofermion occupancy in the corresponding state or subspace, 
the cO and si branches only, and all av branches with finite av pseudofermion occupancy in the corresponding state 
or subspace except the cO and si branches, respectively. 



2. THE PSEUDOFERMION CANONICAL MOMENTUM AND ASSOCIATED FUNCTIONALS 
As mentioned above, the av pseudofermion discrete canonical-momentum values cjj are of the following form, 

5 = qfa) = qj + = ^/r + ; 3 = 1, 2, .... Ku ■ (3) 

where are integers or half-odd integers Q, -/V*„ = N av + N% v , and denotes the number of av pseudofermion 
holes. The latter number equals the corresponding number of av pseudoparticle holes given in Eqs. (B7) and (B8) 
of Ref. 0. Note that besides equaling the number of discrete canonical-momentum values in the av canonical- 
momentum band, N* v = N au + N£ y also equals the number of sites of the av effective lattice which plays an 
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important role in the pseudofermion description. In addition to the av pseudofermions of canonical momentum q, 
there are local av pseudofermions, whose creation and annihilation operators correspond to the sites of the effective 
av lattice. Such a lattice has spatial coordinates Xj = a au j where j — 1, 2, N* v and N*„ is the number of sites 
defined in Eqs. (B.6)-(B.8) and (B.ll) of Ref. Q and aw = L/N*^ is the effective av lattice constant. Each av 
pseudofermion band is associated with an effective av lattice whose length L = N*„ a av is the same as that of the 
original real-space lattice. [|J. The relation between the momentum and local pseudofermion operators is given in 
Eq. (34) of Ref. ®. 

The discrete bare-momentum qj is a good quantum number whose allowed occupancies are one and zero only. 
(Also the corresponding discrete canonical-momentum q~j has allowed occupancies one and zero only.) Thus, for 
the bare-momentum occupancy configuration describing a given energy eigenstate the bare-momentum distribution 
function N av (qj) is such that N av (qj) — 1 for occupied bare- momentum values and N av (qj) = for unoccupied 
bare-momentum values. We denote the ground-state bare-momentum distribution function by N® v (qj). It is given 
in Eqs. (C.1)-(C.3) of Ref. 4]. Although the av pseudoparticles carry bare-momentum qj, one can also label the 
corresponding av pseudofermions by such a bare-momentum. This is because there is a one-to-one correspondence 
between the bare momentum qj and the pseudofermion canonical momentum q~j = qj + Q% v {qj) j L. Thus, when one 
refers to the pseudofermion bare- momentum qj, one means that qj is the bare-momentum value that corresponds 
to the pseudofermion canonical momentum qj = qj + Q® u (qj)/L. The pseudofermion canonical-momentum shift 
functional Q% lJ {qj)/L is given by, 

N* , i 

9%M^Y. E ^, Q ^fe,^)AA QV feO, (4) 



a' if' j 



where, 



AN av (qj) = N ai/ (qj)-N av ( qj ), (5) 



is the av bare-momentum distribution-function deviation. A PS excited energy eigenstate is uniquely defined by 
the values of the set of deviations {AN a „(qj)} for all values of qj corresponding to the av branches with finite 
pseudofermion occupancy in the state and by the values L c _ 1 / 2 and £ S) _i/2- Moreover, the quantity & a v t a'v'(<l,q') 
on the right-hand side of Eq. Q is a function of both the bare-momentum values q and q' given by, 

V« , "'W.?) = *»,<.v jj . jj J) ( 6 ) 

where the function $ a v,a'v'( r i r ') IS the unique solution of the integral equations (A1)-(A13) of Ref. .6). The ground- 
state rapidity functions A9 al/ (q) appearing in Eq. ©, where A® (q) = sinfc°(g) for av = cO, are defined in terms of 
the inverse functions of k°(q) and A°„(g) for v > in Eqs. (A.l) and (A. 2) of Ref. [l5j . 

It is found below that ir ^ au , a'v' (q, s') [° r — T &au, a'v' (?) ?')] i s an elementary two-pseudofermion phase shift such 
that q is the bare-momentum value of a av pseudofermion or av pseudofermion hole scattered by a a'v' pseudofermion 
[or a'v' pseudofermion hole] of bare-momentum q' created under a ground-state - excited-energy-eigenstate transition. 
As discussed in Sec. IV-D, there are no cv ^ cO (and sv ^ si) bare-momentum bands for n — 1 (and m = 0) ground 
states. Indeed, N* v = (and N* v = 0) for such states and then the corresponding ground-state rapidity functions 
A° !/ (g) (and K° sv {q)) cannot be defined. Fortunately, expression © remains valid in that case provided that the 
ground-state rapidity functions are suitably replaced by those of the excited states. (Indeed, we find below that 
the functions 10 are phase shifts originated by well-defined ground-state - excited-state transitions; Thus, in the 
particular case of the n = 1 and/or m — ground states the quantities @l are functionals rather than functions, with 
the rapidity functions for the cv ^ cO and/or sv ^ si branches being those of the excited state under consideration.) 

The form of the functional (QJ reveals that for the initial ground state the discrete canonical-momentum value q~j and 
corresponding discrete bare- momentum value qj are such that q~j — qj . The ground-state continuum bare- momentum 
and canonical- momentum values belong to the domain q € [— <7q„, where the limiting value g°„ reads, 

1% = n ! Qsi = k n ; <l° cv = [tt - 2k F ] , v > ; q° su = [fcp T - k Fl ] , v > 1 . (7) 

The ground-state is described by a compact cO and si pseudofermion finite occupancy for q G [— q Fav , +q F au ], while 
the remaining branches have vanishing ground-state occupancy. Here the cO and si Fermi points are given by, 

IfcO = 2k F 5 q° Fs i = k F i ■ (8) 

Both the limiting values of Eq. and the ground-state Fermi values of Eq. JSJ are given except for corrections of 
order 1 / L. The limiting bare-momentum values and ground-state Fermi bare-momentum including the 1 / L corrections 
are provided in Eqs. (B.14)-(B.17) and (C.4)-(C.ll), respectively, of Ref. Q. 
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3. THE PSEUDOFERMION ENERGY SPECTR UM 



The PS general energy spectrum can be expressed in terms of the pseudofermion canonical-momentum distribution- 
function deviations as follows 0, 

AE = e cQ {q)AM cQ {q) + Y tsi(q) AAfsi(q) 

q=-TT q=-k Fr 

oo +[^ft — ^-Fl] oo +[-7r — 2/cf] 

+ E h + J2 £ e sv (q)AAf s „(q) + E uhb + J2 £ & (q) AN CV (?) . (9) 

f=2 <j=-[/c FT -fc Fl ] w=l q=— [tt— 2/cf] 

Here AJ\f av (q) = AN a „(q) and the energy bands are defined in Eqs. (C.15)-(C.21) of Ref. |4| and plotted in Figs. 
6 to 9 of Ref. for m — 0. The zero-energy level of these energy bands is such that, 

e c0 (±2k F ) = e sl (±k Fl ) = e° cv (±[7t - 2k F ]) = e° sv {±[k F ^ - k n }) = . (10) 

The rotated-electron double occupation D r and the number S r of spin-down rotated-electron singly occupied sites 
whose rotated electrons are not associated with the si pseudofermions play an important role in the finite-energy 
physics and are given by, 

oo oo 

D r = M c , _ 1/2 = [L Ci _ 1/2 + vN cv ] ; S r = [M Sj _ 1/2 - 7V sl ] = [L s ^ _ 1/2 + ^ WV S „] . (11) 

These quantities fully determine the value for the upper-Hubbard band (UHB) energy gap E u hb and spin gap Eh of 
the PS energy spectrum 10 as follows, 

E u hb = 2^ D r ; E h ~2fi HS r . (12) 

Note that E u hb = Eh = for the initial ground state. In equation (|llfl i Qi _iy 2 are the numbers of —1/2 Yang holons 
(a = c) and —1/2 HL spinons {a — s) of the excited energy eigenstate, and N av is the number of av pseudofermions 
of the same state for the av ^ cO, si branches. 



III. THE PSEUDOFERMION SCATTERING THEORY: THE PSEUDOFERMION S MATRIX 

In this section we introduce the pseudofermion scattering theory. We derive the pseudofermion and pseudofermion- 
hole S matrices and discuss their relation to the spectral properties. Our analysis of the problem follows the standard 
quantum non-relativistic scattering theory of spin- less particles |34^ . A ground-state - excited-energy-eigenstate tran- 
sition involves a set of elementary two-pseudofermion scattering events. Such a transition is divided below into three 
steps. The first and second steps have a scatter-less character and lead to the excitation momentum and energy asso- 
ciated with the transition. Through these two steps the ground-state - excited-energy-eigenstate transition reaches the 
many-pseudofermion "in" state, which contains the one-pseudofcrmion "in" asymptote states of the pseudofermion 
scattering theory. The third step corresponds to a well defined set of elementary two-pseudofermion scattering events 
which give rise to the "out" state and conserve both the momentum and the energy. The latter many-pseudofermion 
state is a excited energy eigenstate and contains the one-pseudofermion "out" asymptote states of the pseudofermion 
scattering theory. An important point for applications to the study of the finite-energy spectral and dynamical 
properties is that all "in" and "out" states of the theory are energy eigenstates 

We start our analysis by a discussion of the PS subspaces and the pseudofermion creation and annihilation operator 
anti-commutation relations. 



A. PS SUBSPACES AND THE PSEUDOFERMION OPERATOR ANTICOMMUTATORS 



Several PS subspaces play an important role in the pseudofermion theory. An electronic ensemble space is a subspace 
spanned by all energy eigenstates with the same values for the electronic numbers ATj and N±. A CPHS ensemble 
space is a subspace spanned by all energy eigenstates with the same values for the numbers {M a ±i/ 2 } of ±1/2 holons 
(a = c) and ±1/2 spinons (a = c) (In CPHS ensemble space, CPHS refers to cO pseudofermion, holon, and 
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spinon.) A CPHS ensemble subspace is spanned by all energy eigenstates with the same values for the sets of numbers 
N c0 , {N cu }, {N SI/ }, L c> _i/ 2 , and L s ,_i/ 2 such that v= 1,2, .... 

At zero absolute temperature the pseudofermion description corresponds to a ground-state normal-ordered theory 
@ . Thus, there is a pseudofermion theory for each initial ground state. The minimum excitation energy value of the 
energy eigenstates that span a given CPHS ensemble subspace involves the gap parameters of Eq. (jT2"|) and is given 
by 

lu q = LU Q (D ri S r ) = E uhb +E h = 2(i D r + 2/io H S r . (13) 

For the ground state D r = S r = and thus ujo = 0. The application onto the latter state of an one-, two-, or any 
other finite-number-electron operator O' generates an excitation which can be described as a suitable superposition 
of PS excited e nerg y ei gens tates. The PDT provides the matrix elements associated with the coefficients of such 
a superposition [l4 ITM Tl6|| . Therefore, the elementary ground-state - excited-energy-eigenstate transition plays a 
central role in the pseudofermion theory 

An excitation 0'\GS) associated with small values for the deviations ATV-j- and AN^ is contained in a well defined 
direct sum of CPHS ensemble subspaces, 

^cphs © Scphs © Scphs © Scphs © — ■ (14) 

Here S l cphs with i = 1, 2, 3, ... corresponds to different CPHS ensemble subspaces. The pseudofermion, Yang holon, 
and HL spinon number deviations of all the CPHS ensemble subspaces of such a direct sum obey the sum rules (18) 
and (19) of Ref. 

and the selection rules given in Eq. (21) of the same reference. 
Let us consider a av pseudofermion of canonical momentum q and a aV pseudofermion of canonical momentum q' 
such that the canonical-momentum values q and q' correspond to an excited energy eigenstate and the initial ground 
state, respectively. Following the results of Refs. [Til ITEl Il6[. the anticommutators involving the creation and/or 
annihilation operators of these two pseudofermions play a key role in the study of the finite-energy spectral and 
dynamical properties. Such anticommutators read |(J, 

s sinf Q au (q)/2 ) 

I ft f-, I _ °av,a'v' j(g-g')/2 iQ a „M/2 \ /. r f t f t \ - I f- f-, , ,\ - ft 

XJq, avi Jq',a'v'J — ^ e e sm([q — q']/2) ' q ' QI " a ' v ' ~ aly ' °" / J — u • V °> 

Here Q a u(q)/^ is the value of the following functional for the above excited energy eigenstate, 

<M<&)/2 = Q° av /2 + Qt(q,)/2 , (16) 

where Q® u (qj)/2 is the functional given in Eq. (J1J whose bare-momentum distribution function deviations correspond 
to that state and 0, 

oo oc 

Q°co = I] J2 ANa » even; Q°o = ±^; Y, ANa » odd ; 

OL — C. S V — l Oi — C. s V — l 

Ql u = 0; AN c0 + AN a „ even; Q° av = ±n ; AN c0 + AN a „ odd; a = c, s , v>ft. (17) 

When Q® al , = ±7r for the av ^ cO bands, the uniquely chosen and only permitted value Q av — vr or Q\ v = —it 
is that which leads to symmetrical limiting discrete bare-momentum values ±[ir/L] [N* v — 1] for the excited-state 
bare-momentum band. (See Eq. (B.14) of Ref. Q.) In turn, for the cO branch the bare-momentum band width is 
2ir. Thus, in this case = n an d Qco = — w l ea d to allowed occupancy configurations of alternative excited energy 
eigenstates. (In the particular case that the cO band is full for the excited energy eigenstate, the two values = n 
and Q° c0 = —ir refer to two equivalent representations of that state.) 

The quantity Q° av /L is the shift in the discrete bare-momentum value qj = [27r/i]J" ly of Eq. © that arises due to 
the transition from the ground state to the excited energy eigenstate. Furthermore, Q a ^(qj)/L is the corresponding 
shift in the discrete canonical-momentum values that occurs as a result of the same transition. It follows from Eq. 
(I15|l that the functional Q a v{q)/2 fully controls the pseudofermion anticommutators associated with the ground-state 
- excited-energy-eigenstate transition. The one- and two-electron matrix elements between the initial ground state 
and the excited energy eigenstates can be expressed in terms of the anticommutators l|15|) |14| . This justifies the 
importance of the functional Q av {q)/2 given in Eq. 11611 . once it controls the quantum overlaps associated with the 
one- and two-electron finite-energy spectral properties |l4 ITU IT"6| . 

Let us introduce the av pseudofermion scattering theory. Within such a theory the functional Q a v(q)/2 is an 
overall av pseudofermion or hole phase shift. 
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B. THE GROUND-STATE - VIRTUAL-STATE TRANSITION 



From now on and until section IV-D our analysis refers to initial ground states with density values in the ranges 
< n < 1 and < n < m. The specific properties of the scattering theory for initial ground states corresponding to 
n = 1 and/or m — are considered in that section. The preliminary analysis of the pseudofermion scattering problem 
presented in Ref. [2!| divided each transition from the initial ground state to a PS excited energy eigenstate into two 
main steps. However, it is useful for our study to divide the first step considered in that reference into two processes. 
The first process considered here is a scatter-less finite-energy and finite-momentum excitation which transforms the 
ground state onto a well defined virtual state. For v > branches, that excitation can involve a change in the number 
of discrete bare-momentum values given by, 

oc oc 

AN* t = AN c0 - AN sl -2^2 AN s» \ AiV*„ = AL a + 2 (z/ - v)AN av , ; av £ cO, si . (18) 

v=2 v'=v+l 

For the initial ground state these numbers read, 

N° w * = (N a - N) ; N s {*=N r> = (iV T - N L ) , v>l, (19) 

and N®q* = N* is given by N* = N a for the whole Hilbert space. Although the av ^ cO, si branches have no finite 
pseudofermion occupancy in the initial ground state, for densities < n < 1 and < m < n one can define the values 
N*„ — N^ v for the corresponding empty bands. For the av ^ cO, 1, si branches, those are the numbers N®f and 
N^* given in Eq. 1)19(1 . Thus, for av ^ cO, si branches with finite pseudofermion occupancy in the virtual state the 
deviations i|18|) and discrete bare-momentum shifts l|17fl are relative to the values of these empty bands. 

In addition and following the change in the number of discrete bare-momentum values, this excitation also involves 
the pseudofermion creation and annihilation processes and pseudofermion particle-hole processes associated with PS 
excited states. The momentum and energy of this ground-state - virtual-state transition is given by, 

AP = ^N c0 (q)+ q^sl(q)+Yl E 1 A NsA<l) 

q— — 77 q— — k F ^ v—2 q— — [fc FT —^fi] 

oo oo +[tt— 2kp] 

+ * [La, -1/2 + E VN ^\ + £ Yj ^ ~ «] AiV -(«) > ( 2 °) 
v=l v=\ q= — [tt— 2fc F ] 

and 

AE = eco(g)AJVco(g) + £ e sl (q) AN sl (q) 

q=-TT g=-fc FT 

oo +[kFy-k Fl ] oo +[-!r-2k F ] 

+ E h + J2 E e° s „(<l)&N sl/ (q)+E uhb +J2 £ e° c „(q) AN c „(q) , (21) 

v—2 q=-[k F} -fcjrj.] »=1 q=-[ir-2k F ] 

respectively, where all quantities were defined above. Except for 1/L energy corrections, the energy spectra (0 and 
(|21|l are identical. In this first scatter-less step the pseudofermions acquire the excitation momentum and energy 
needed for the second and third steps. 



C. PSEUDOFERMION SCATTERING PROCESSES, S MATRICES, AND PHASE SHIFTS 

In order to study the second and third processes of the ground-state - excited-energy-eigenstate transition, it is 
useful to express the many-pseudofermion states and operators in terms of one-pseudofermion states and operators, 
respectively. The PS energy and momentum eigenstates can be written as direct products of states spanned by 
the occupancy configurations of each of the av branches with finite pseudofermion occupancy in the state under 
consideration. Moreover, the many-pseudofermion states spanned by occupancy configurations of each av branch 
can be expressed as a direct product of N*„ one-pseudofermion states, each referring to one discrete bare-momentum 
value q jt where j = 1, 2, ...,N* V . 
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Within the pseudofermion description, the ID Hubbard model in normal order relative to the initial ground state 
has no residual-interaction energy terms. Thus, when acting in the PS it has a uniquely defined expression of the 
general form Q, 

av j—1 a 

where we denoted the ground-state normal ordered Hamiltonian by : H :, H avq . is the one-pseudofermion Hamiltonian 
which describes the av pseudofermion or hole of bare-momentum qj , and H a refers to the Yang holons (a — c) and 
HL spinons (a — s), which are scatter-less objects. Thus, for each many-pseudofermion PS virtual state reached in 
the first step of the transition from the ground state to the excited energy eigenstate, the number of Hamiltonians 
H av ^ qj equals that of one-pseudofermion states of the virtual state given by, 

K Q + N* sl + 0(\AN av \)N* v , (23) 

ai/^cO, si 



(XV 



Here 6{x) = 1 for x > and 9(x) = for x = and N* = N c0 + N% , N* x = N sl + N> s \, and N* u = N av + N ( 
refer to the virtual state and corresponding excited energy eigenstate under consideration. The pseudofermion-hole 
numbers of these expressions read 

oo oo 

N^ = N a -N c0 ; Nl = N%-Y i (y + i/-\v-^N a ,', N* v = iV c0 - ^ (y + v' - \v - v'\)N sv , . (24) 

u' = l v' = \ 

For the pseudofermion description only momentum and energy contributions of order zero and one in 1/L are 
physical and thus our analysis refers to periodic boundary conditions and large values of L such that L » 1. The 
second scatter-less process generates the "in" state. Indeed, the one-pseudofermion states belonging to the many- 
pseudofermion "in" state are the "in" asymptote states of the pseudofermion scattering theory. The generator of the 
virtual-state - "in" -state transition is of the form, 

s°=nn . ( 2s ) 

QLV j — 1 

where S^ v q . is a well-defined one-pseudofermion unitary operator. Application of S^ v q . onto the corresponding 
one-pseudofermion state of the many-pseudofermion virtual state shifts its discrete bare-momentum value qj to the 
bare-momentum value qj + Q% v /L, where Q° av is given in Eq. I|17() . 

Finally, the third step consists of a set of two-pseudofermion scattering events. It corresponds to the "in" -state - 
"out"-state transition, where the latter state is the PS excited energy eigenstate under consideration. The generator 
of that transition is the following operator, 

^ = n n § t, Qi > (26) 

ay j — 1 

where S% v is a well-defined one-pseudofermion scattering unitary operator. The one-pseudofermion states belonging 
to the many-pseudofermion "out" state are the "out" asymptote pseudofermion scattering states. Application of 
Stv,qj 0lrt0 the corresponding one-pseudofermion state of the many-pseudofermion "in" state shifts its discrete bare- 
momentum value qj + Q° au / L to the "out"-state discrete canonical-momentum value qj + Q a v{<lj) / L. It follows that 
the generator of the virtual-state - "out" -state transition is the unitary operator, 

s = s*s° = n n - ( 2? ) 

QLU j — 1 

where S av ^ qj is the one-pseudofermion or hole unitary S a ^^ qj — S^ vq .S% vq . operator. Application of the latter 
operator onto the corresponding one-pseudofermion state of the many-pseudofermion virtual state shifts its discrete 
bare-momentum value qj directly to the "out"-state discrete canonical-momentum value qj + Q a v(.q_j) I L. 



10 



The virtual state, "in" state, and "out" state are PS excited energy eigenstates, as further discussed below. Thus, 
that the one-pseudofermion states of the many-pseudofermion "in" state and "out" state are the "in" and "out" 
asymptote pseudofermion scattering states, respectively, implies that the one-pseudofermion Hamiltonian H av qj plays 

the role of the unperturbed Hamiltonian Hq of the spin-less one-particle nonrelativistic scattering theory j^ij . Indeed, 
the unitary S av ^ qj operator (and the scattering unitary S^ v q . operator) commutes with the Hamiltonian H av ^ qj and 
thus the one-pseudofermion "in" and "out" asymptote scattering states are energy eigenstates of H av ^ qj and eigenstates 

of Saurfj (and S^ uq .). It follows that the matrix elements between one-pseudofermion states of S ail , qj (and S^ vq .) 
are diagonal and thus these operators are fully defined by the set of their eigenvalues belonging to these states. The 
same applies to the generator S (and S^) given in Eq. i(27() (and Eq. 1(26(1). The matrix elements of that generator 
between many-pseudofermion virtual states (and "in" states) are also diagonal and thus it is fully defined by the set 
of its eigenvalues belonging to such states. Importantly, the virtual state and the "in" state of a given ground-state 
transition correspond to the same excited energy eigenstate of the ID Hubbard model, which by construction is the 
"out" state. Indeed, the virtual state, "in" state, and "out" state only differ by mere overall phase factors whose 
general functional expression is given below. That the many-pseudofermion "in" and "out" states which are a direct 
product of one-pseudofermion "in" and "out" asymptote pseudofermion scattering states, respectively, are PS excited 
energy eigenstates plays a major role in the pseudofermion scattering theory. 

Since S^ q . and S av ^ qj are unitary, each of their eigenvalues has modulus one and can be written as the exponent 
of a purely imaginary number given by, 

N*, u , 

StM = e i0 -fe) = ]J f[ S w , aV fe (9j ,); j = l,2,..,Ar* v 

ol'v' j' — l 

S av { qj ) = e ^to) =e «3° wS *(^) ; j = l,2,...,N* v . (28) 

Here Qavilj) an d Qav{qj) are the functionals defined by Eqs. Q and (|16|l . respectively. By use of the former 
functional we find that, 

S„,oyfe,^) = ^»«,^(*.V)^M , (29) 

where the functions 7r <$> av , a 'v' {Qj, Qj') are given in Eq. JBJ. The main point is that except for the occupancy 
configuration changes produced by the ground-state - virtual-state transition, the only effect of under a ground-state - 
excited-energy-eigenstate transition moving the av pseudofermion or hole of initial virtual-state canonical-momentum 
qj = qj once around the length L lattice ring is that its wave function acquires the overall phase factor S av (qj) given 
in Eq. 1)28(1 . This property follows from the form of the energy spectrum of the pseudofermions, which in contrast to 
that of the corresponding pseudoparticles of Ref . Q has no residual interaction terms . 

The phase factor S aI/ta > U '(qj,qj>) of Eq. (|29|l in the wave function of the av pseudofermion or hole of bare- 
momentum qj results from an elementary two-pseudofermion zero-momentum forward-scattering event whose scat- 
tering center is a aV pseudofermion (AA^v/ (qy) = 1) or a'v 1 pseudofermion hole (AJV Q! / I/ '(g,v) = — 1) created 
under the ground-state - excited-state transition. Thus, the third step of that transition involves a well-defined 
set of elementary two-pseudofermion scattering events where all av pseudofermions and av pseudofermion holes of 
bare-momentum qj + Q a av /L of the "in" state are the scatterers, which leads to the overall scattering phase factor 
S®v(qj) in their wave function given in Eq. ((28(1 . That the scattering centers are the a'v 1 pseudofermions or pseud- 
ofermion holes of bare momentum qj' + Q av l L created under the ground-state - "in" -state transition is confirmed by 
noting that S al/ ^ a ' v i(qj,qji) = 1 for AN a r„r(qji) — 0. Thus, out of the scatterers whose number equals that of the 
one-pseudofermion states given in Eq. 1(23(1 , the scattering centers are only those whose bare- momentum distribution- 
function deviation is finite. However, out of the above scatterers, only a subclass of scatterers contributes significantly 
to the spectral properties 0, . 

The following properties play an important role in the pseudofermion scattering theory: 

1. The elementary two-pseudofermion scattering processes associated with the phase factors 1(29(1 conserve the total 
energy and total momentum. Such an energy conservation is further discussed in Appendix A. 

2. The elementary two-pseudofermion scattering processes are of zero-momentum forward-scattering type and thus 
conserve the individual "in" asymptote av pseudofermion momentum value qj + Q® al/ / L and energy. 

3. These processes also conserve the av branch, usually called channel in the scattering language |34j . 

4. The scattering amplitude does not connect quantum objects with different r\ spin or spin. 
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5. For each av pseudofermion or pseudofermion hole of virtual-state bare-momentum qj, the S matrix associated 
with the ground-state - excited-energy-eigenstate transition is simply the phase factor S au (qj) given in Eq. I|28() . 

The one-particle phase factor si(E) of Eq. (6.8) of Ref. [34[ whose expression is given in Eq. (6.9) of the 
same reference corresponds to the one-cw-pseudofermion or hole phase factor S^ [/ (qj) with the energy E and the 
quantum numbers / and m replaced by the bare- momentum qj. Indeed, while the av pseudofermion or hole energy 
is uniquely defined by the absolute value \qj\, in ID the sign of qj corresponds to the three-dimensional angular- 
momentum quantum numbers. Another difference is that si(E) is associated with a single scattering event whereas 
Sauili) results in general from several scattering events. Each of such events corresponds to a well defined factor 
S a v, a'v' {UjjQj' ) °f form 1 (2*5(1 in the expression of S® v (qj) given in Eq. 1(28)1. There are as many of such factors as aV 
pseudofermion and hole scattering centers created under the transition to the virtual state and corresponding excited 
energy eigenstate under consideration. The factor 2 in the phase factor of Eq. (6.9) of Ref. [34( corresponds to the 
phase-shift definition of the standard nonrelativistic scattering theory for spin-less particles. As discussed below, we 
use here such a definition which introduces the overall scattering phase shift S® v (qj) — Qauilj)/^ an d overall phase 
shift 5 av (q,j) — Q a y{qj) /2. However, if instead we insert a factor 1, we would introduce the overall scattering phase 
shift QuviQj) an d overall phase shift Q av (qj) whose values are defined only to within addition of an arbitrary multiple 
of 2n. That is the choice for the phase shift definition used in Refs. [3ll l32L l33| . For our definition the phase shifts 
are instead given only to within addition of an arbitrary multiple of tt, as in Ref. |34| . 

Application of the unitary S av ^ qj operator onto its one-pseudofermion state of the many-pseudofermion virtual 
state generates the corresponding one-pseudofermion state of the many-pseudofermion "out" state. The latter one- 
pseudofermion state equals the former one multiplied by the phase factor S au (qj) of Eq. (|28[l . (Applying the scattering 
unitary 5* operator onto its one-pseudofermion state of the many-pseudofermion "in" state also generates the 
corresponding one-pseudofermion state of the many-pseudofermion "out" state; The latter one-pseudofermion state 
equals the former one multiplied by the phase factor S®„(qj) of Eq. J2SJ-) It follows that the many-pseudofermion 
virtual states (and "in" states) are eigenstates of the generator S (and S^) given in Eq. (|27|) (and Eq. The 
eigenvalue of S belonging to a PS virtual state is given by, 

N*„ AT*„ 

St = ^ 25t =Y[f[S av (q j ); £ T = ^]TQa„fe)/2. (30) 

av j — 1 olu j—1 

Thus, the corresponding "out" state equals the virtual state multiplied by the phase factor St- Furthermore, the 
eigenvalue of S* belonging to a PS "in" state is given by, 

N*„ N* u 

5* = e i2S * = [] I! S t(<li) 5 *?=EE • (31) 

av j — 1 av j — 1 

Again, the corresponding "out" state equals the "in" state multiplied by the phase factor 5*. Since the "out" state 
is by construction an energy eigenstate of the ID Hubbard model, this result confirms that the corresponding virtual 
and "in" states are also energy eigenstates of the model: the "out" state only differs from the latter two states by 
the phase factors given in Eqs. H3U[) and (|31fl . respectively. The general expressions """J and (|16|l for the functionals 
Qttvilj) an d Qav{qj) define uniquely the eigenvalues 5 1 * and St of S^ and S for the whole set of "in" states and 
virtual states, respectively, corresponding to the excited energy eigenstates that span the PS. 

The factorization of the BA bare S matrix for the original spin 1/2 electrons is associated with the so called Yang- 
Baxter Equation (YBE) [30l l32|. On the other hand, the factorization of the av pseudofermion or hole S matrix 
S au {qj) given in Eq. 1)28(1 . in terms of the elementary two-pseudofermion S matrices Sau,a'u'(QjtQj')j Eq. (|29|l . is 
commutative. Such commutativity is stronger than the symmetry associated with the YBE. The pseudofermion S 
matrix commutative factorization is required by the form of the pseudofermion occupancy configurations that describe 
the PS excited energy eigenstates. These states are direct products of the one-pseudofermion scattering states of the 
theory and are described by well defined occupancy configurations of rotated electrons. All such configurations have the 
same number of rotated-electron occupied sites, unoccupied sites, spin-up singly occupied sites, and spin-down singly 
occupied sites. However, the relative position of these quantum objects is different in each occupancy configuration. 
There is a one-to-one correspondence between these rotated-electron configurations and the local pseudofermion, 
Yang holon, and HL spinon occupancy configurations that describe the same state. Again, the number of local av 
pseudofermions belonging to av branches with finite occupancy in the virtual state are the same for all occupancy 
configurations. As for the rotated electrons, the relative position of these quantum objects is different in each 
configuration. Thus, when under a specific ground-state - excited-energy-eigenstate transition a av pseudofermion 
or hole moves around the lattice ring, it scatters the same a'v' pseudofermion or hole scattering centers, but in 
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different order for different occupancy configurations. However, it is required that the phase factor e^""™' acquired 
by the av pseudofermion or hole should be the same, independently of that order. This implies the commutativity 
of the S- matrix factors S a v,a'v'{<lji<lj') m the overall S matrix S av (qj) whose expression is given in Eq. J2HJ- Such 
commutativity follows from the elementary S matrices S a i>,a'v'(qj,qj') given in Eq. Ij29(l being simple phase factors, 
instead of matrices of dimension larger than one. This seems to be inconsistent with all energy eigenstates being 
described by occupancy configurations which, besides cO pseudofermions, involve finite-spin 1/2 spinons and ?7-spin 
1/2 holons Indeed, the 5* matrix of finite- 77-spin or spin quantum objects is a matrix of dimension larger than one 
[33ll34l| . However, in spite of these fmite-77-spin and finite-spin objects, due to symmetry the system self organizes in 
such a way that, in addition to the 77-spin-less and spin-less cO pseudofermions, the scatterers and scattering centers 
are the 77-spin zero 2^-holon composite cv pseudofermions and spin zero 2i/-spinon composite sv pseudofermions. 
Symmetry requirements also imply that the 77-spin 1/2 Yang holons and spin 1/2 HL spinons 0,0 are purely scatter- 
less objects, as discussed in the following. 

Thus, it is the 77-spin-neutral (and spin-neutral) character of the 2i/-holon (and 2i/-spinon) composite pseudofermion 
scatterers and scattering centers and the 77-spin-less and spin-less character of the cO pseudofermion scatterers and 
scattering centers which is behind the dimension of their S matrix S au (qj) given in Eq. (J2SJ. Interestingly, in the 
pseudofermion scattering theory the relation of the composite cv pseudofermion (and sv pseudofermion) scatterers 
and scattering centers to the holons (and spinons) has similarities with that of the composite physical particles to 
the quarks in chromo dynamics [35^ . Within the latter theory the quarks have color but all quark-composite physical 
particles are color- neutral. Here the holons (and spinons) have 77 spin 1/2 (and spin 1/2) but the 2^-holon (and 
2^-spinon) composite pseudofermion scatterers and scattering centers are 77-spin-neutral (and spin- neutral) . 

In turn, the holon-holon and spinon-spinon S matrices of the conventional spinon-holon scattering theory of Refs. 
[30t l3ll l32| do not have the above commutative properties. Indeed, within that theory the scatterers and scattering 
centers have ?7-spin 1/2 or spin 1/2 and for initial ground states with densities in the ranges < n < 1 and < m < n 
many one-particle scattering states do not refer to energy eigenstates [33j. (The above mentioned requirement for 
commutative factorization of the S matrix applies when the one-particle scattering states belong to many-particle 
energy eigenstates and the scatterers and scattering centers are 77-spin-neutral and/ or spin-neutral.) On the other 
hand, for the n = 1 and m = initial ground state considered in Refs. [3(], the one-particle scattering 

states of the conventional spinon-holon representation of these references belong to many-particle energy eigenstates 
33]. However, since the scatterers and scattering centers of that theory have 77-spin 1/2 or spin 1/2, instead of the 
commutative factorization it is required that the S matrix has a YBE like factorization, as the BA bare S matrix of 
the original spin 1/2 electrons. 



D. SCATTERERS, SCATTERING CENTERS, AND SYMMETRY 

The three generators of both the 77-spin and spin SU (2) algebras commute with the electron - rotated-electron 
unitary operator 0,0. As a result, the symmetry transformations of the av pseudofermions, ±1/2 Yang holons, and 
±1/2 HL spinons play an important role in their scattering properties. 

For initial ground states with densities such that < n < 1 and < m < n there is only finite occupancy 
for cO pseudofermions and holes, si pseudofermions and holes, ±1/2 Yang holons, and ±1/2 HL spinons. The 
corresponding occupation numbers read Y c0 = N, Y' l = [N a - N], N sl = iVj,, = [iVf - iVjJ, £ c ,+i/2 = [N a - N], 
and L S) +i/2 = [Nf — iVjJ. However, that for such ground states N£ = L c ,+i/2 and iVj\ = £ s ,+i/2 does not imply 
that the cO pseudofermion holes and the si pseudofermion holes are the same quantum objects as the ±1/2 Yang 
holons and ±1/2 HL spinons, respectively. For instance, the [N a — N] ±1/2 Yang holons have finite 77 spin 1/2 and 
the [iVf — YJ ±1/2 HL spinons finite spin 1/2 and both these objects are dispersion-less and have zero momentum 
and energy. In contrast, the cO pseudofermion holes are 77-spin-less and spin-less and the si pseudofermion holes 
have spin zero and are ?7-spin-less and both these types of pseudofermion holes have a momentum-dependent energy 
dispersion, e c o(q) and e s i(q), respectively. The ground-state canonical-momentum distributions of the [N a — N] cO 
pseudofermion holes and [Yj — YjJ si pseudofermion holes correspond to compact domains such that 2Uf < \q\ < tt 
and kpi < \q'\ < kp^, respectively. Moreover, while under a ground-state - excited-energy-eigenstate transition the 
cO and si pseudofermion-hole discrete canonical- momentum values q~j and q'j acquire a shift given by Q c o(qj)/ L and 
Q s i(q' j)/L, respectively, the momentum zero of the ±1/2 Yang holons and ±1/2 HL spinons remains unchanged. 
Also the momentum values 7r and zero of the —1/2 Yang holons and —1/2 HL spinons, respectively, remain unchanged 
under such a transition. 

The form of the scattering part of the overall phase shift I jltjl . Eq. (0}, reveals that the value of such a phase-shift 
functional is independent of the changes in the occupation numbers of the ±1/2 Yang holons and ±1/2 HL spinons. 
Thus, these objects are not scattering centers. Moreover, they are not scatterers, once their momentum values remain 
unchanged under the transition from the "in" state to the "out" state (excited energy eigenstate). Such a scatter- less 
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character of the ±1/2 Yang holons and ±1/2 HL spinons is related to the above symmetry. Indeed, the ±1/2 Yang 
holons (and ±1/2 HL spinons) are created and annihilated by the 77-spin (and spin) off-diagonal generators. Since 
these generators commute with the electron - rotated-electron unitary operator, they have the same [/-independent 
expressions both in terms of electronic and rotated-electronic operators 0,0]. That the ±1/2 Yang holons and ±1/2 
HL spinons are neither scatterers nor scattering centers is consistent with that property. 

In contrast, since the pseudofermions and pseudofermion holes are not in general invariant under the electron - 
rotated-electron unitary transformation jj, |^| , their creation and annihilation operators have U -dependent expressions 
in terms of electronic operators. That these quantum objects are scatterers and scattering centers is consistent with 
such a property. Moreover, these objects are not transformed by the 77-spin and spin generators. This latter symmetry 
is behind the cO pseudofermion being a 7/-spin-less and spin-less object and for v > the av pseudofermions being 
77-spin (a = c) and spin (a — s) singlet 2^-holon and 2i/-spinon composite objects, respectively. 

The particular case of the invariance under the electron - rotated-electron unitary transformation of av 7^ cO, si 
pseudofermions of limiting canonical-momentum values q = ±g°„ and the corresponding separation of these objects 
into 2v independent holons {av = cv) or spinons (av = sv) and av 7^ cO, si FP (Fermi-po'mt) current scattering 
centers is discussed below in Sec. IV-C. 



E. PSEUDOFERMION 5 MATRIX AND THE FINITE-ENERGY SPECTRAL PROPERTIES 



The simple form of the pseudofermion and hole S matrix renders the pseudofermion description particularly suitable 
for the study of the unusual finite-energy spectral properties of the model. Such a simplicity of the S matrix form 
results in part from all "in" and "out" states being energy eigenstates. Indeed, the latter states can be expressed 
as direct products of the "in" and "out" one-pseudofcrmion scattering states of the theory. Fortunately, this allows 
the use of Lehmann representations for the study of the spectral functions |14lll5lll6l |. Moreover, the simple form of 
the pseudofermion and hole S matrix simplifies the evaluation of the spectral-function matrix elements between the 
initial ground state and the exact excited energy eigenstates, which are the "out" states. Indeed, the anticommutation 
relations (|15|) can be expressed solely in terms of the difference q — q' of the two pseudofermion momenta and the S 
matrix S al ,(q) given in Eq. (|28|l of the pseudofermion associated with the excited energy eigenstate as follows, 



{fq,avi fq',a'v'} — S av (c[) €. ^ q ^ . , (32) 



JV* 



1/2 , ^ 

-i{q-q')/2 



S a „(q) 



1/2 



«?-9']/2) 



and the anticommutators between two av pseudofermion creation or annihilation operators vanish. This reveals that 
the S matrix S aiy {qj) fully controls the pseudofermion anticommutators. Since within the PDT these anticommutators 
determine the value of the matrix elements between the ground state and the excited energy eigenstates |l4l Il5t 1 1 6| , 
it follows that the S matrix S av (qj) controls the finite-energy spectral properties. If it had dimension larger than 
one, the pseudofermion algebra would be much more involved, for the pseudofermion anticommutators would also 
be matrices of dimension larger than one. The problem of the evaluation of the spectral-function matrix elements 
between energy eigenstates simplifies for the pseudofermion representation because the PS subspaces associated with 
a given one- or two-electron spectral function can be expressed in terms of direct products corresponding to each of the 
av branch quantum- number occupancy configurations of branches with finite pseudofermion occupancy 0, ITEl Il6| . 
For these matrix elements the direct product is associated with the commutative factorization of the S matrix l|28|) 
in terms of the elementary S matrices S aUya i v >(qj,qji), Eq. (|29(l . Therefore, the use of the pseudofermion description 
considerably simplifies the study of the exotic quantum- liquid finite-energy spectral properties 0, 0, |23l 0, H3] • 

The PDT of Refs. [TH confirms that the overall av pseudofermion and hole phase-shift functional ll(i|) 
associated with the S matrix S au (qj) of Eq. (|28|l fully controls the one- and two-electron spectral properties through 
the pseudofermion anticommutators. An one- or two-electron excitation 0^\GS) is contained in a well defined direct 
sum (|14fl of CPHS ensemble subspaces and is a superposition of excited energy eigenstates. Each of these states is 
described by a set of deviations \ AN nu (q^ ) }. Eq. JSJ), involving a finite number of scattering centers. Such deviations 
obey the sum rules (18) and (19) of Ref. [14j and the selection rules given in Eq. (21) of the same reference. 

One can consider that the order of the direct sum (|14|) of CPHS ensemble subspaces associated with a given one- 
or two-electron excitation is that of the increasing values of the minimum energy w {D ri S r ), Eq. lfH?)l . of each CPHS 
ensemble subspace relative to the initial ground state. For fixed value of S r and D r , these energy values are such 
that, 

Wo(0, S r ) < wo(l, S r ) < L0 (2, S r ) < ... , (33) 

and 



tu {D r , 0) < u> (D r , 1) < CJoPr, 2) < ... , 



(34) 
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respectively. Note that for fixed values of electronic density n and spin density m, the ordering of the energies 
Wo(D r , S r ) is well defined and such that, 

^ (0, 0) < min{ Wo (0, 1), w Q (l, 0)} < ... . (35) 

An important property is that for a given energy range < AE < loq the number of CPHS ensemble subspaces 
of the direct sum 1|14|) is finite. Moreover, as each energy u>o(D r , S r ) is reached, an increasing number of channels 
open up that correspond to CPHS ensemble subspaces of increasing energy. By new channels we mean here av 
pseudofermion branches of increasing v value. Thus, for a given energy range < AE < luq there is a well defined set 
of pseudofermion S matrices, whose "in" and "out" asymptote one-pseudofermion scattering states correspond to the 
set of excited energy eigenstates spanning the direct sum (I14fl of CPHS ensemble subspaces of energy smaller than or 
equal to u>o- The pseudofermion S matrices belonging to that finite set control the spectral properties for energies in 
the above range < AE < ujq. 

The pseudofermion representation is valid for L >> 1. Indeed, the expressions of all quantities derived by use of that 
representation are physical up to first order in 1/L. Thus, for the study of some properties one can replace the discrete 
bare-momentum qj and canonical-momentum q~j by a continuum bare-momentum variable q and canonical-momentum 
variable q, respectively. The ground-state is then described by a cO occupancy for \q\ < 2kp and unoccupancy for 
2kp < \q\ < it and a si occupancy for |g| < kpi and unoccupancy for kpi < \q\ < kpp (All remaining av 
pseudofermion bands are unoccupied for the ground state.) Thus, in the continuum momentum limit, the global 
canonical-momentum shift Q av {q)/L — Q Q av /L + Q* J/ (g)/L contributes to the spectral properties mainly through 
the av = cO, si branches for q values in the vicinity of the Fermi points iq^n, as confirmed by the studies of 
Refs. GHHHl. However, often such a limit must be taken in the end of the calculations. Otherwise, one would 
loose the information contained in the overall pseudofermion or hole phase shifts studied below, which correspond to 
canonical- momentum shifts, Q a „(q)/L, of the order of 1/L. 



IV. THE PSEUDOFERMION SCATTERING THEORY: PSEUDOFERMION PHASE SHIFTS 

In this section we study the av pseudofermion phase shifts associated with the S matrix introduced above. 



A. PHASE-SHIFT DEFINITION 

The effective av lattices have the same length L as the rotated-electron and electronic lattice [jj. As above, our 
analysis refers to periodic boundary conditions and L » 1. Depending on the asymptote coordinate choice, there are 
two possible definitions for the av pseudofermion phase shifts associated with the S matrix S au {qj) given in Eq. J2HJ>- 
The choice of either definition is a matter of taste and the uniquely defined quantity is the S matrix. The two choices 
of asymptote coordinates for the av pseudofermion or av pseudofermion hole correspond to x G (—L/2, +L/2) and 
x £ (0, +L). 

If when moving around the lattice ring the av pseudofermion (or hole) departures from the point x = —L/2 and 
arrives to x = L/2, one finds that, 

lim qx = qx + Q° a j2 + Q* v (q)/2 = qx + 5 av (q), (36) 

x— *L/2 

where 

S m ,(q) = Q m ,(q)/2. (37) 

For this asymptote coordinate choice, 8 au {q) is the overall av pseudofermion or hole phase shift whose value is given 
only to within addition of an arbitrary multiple of tt. Moreover, from analysis of Eqs. and (|16fl it follows that 
'R'&av,a'v'(qj,qj') is an elementary two-pseudofermion phase shift. This phase-shift definition corresponds to that 
used in standard quantum non-relativistic scattering theory |34|. such that the S matrix S au (qj) given in Eq. I|28|) 
can be written as, 

S av ( qj ) = e i25 °">^ , j = 1, 2, N* av . (38) 

The factor 2 appearing in the exponential argument of Eq. (|38f) corresponds to the usual form of the S matrix for 
that theory. In reference j3|| it is found that such a phase-shift definition is consistent with an exact theorem due to 
Fumi 37]. Note that for the phase-shift definition i|37|). Q%„/2 — 0, =F 7r /2 corresponds to the scatter-less term —lir/2 
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of the three-dimensional partial- wave problem of angular momentum I |3J, |3?| ■ Although the angular momentum 
vanishes in ID, the scatter-less phase shift <|17|) plays a similar role. Here we follow the definition of the standard 
quantum non-relativistic scattering theory and choose the overall av pseudofermion phase shift definition Q av (q)/2 
associated with Eq. {35} . 

In turn, if when moving around the lattice ring the av pseudofermion (or hole) departures from the point x = 
and arrives to x = L, one finds that, 

lim qx = qx + Q° av + Q* v {q) = + Q av (q) , (39) 

x — >L 

where q refers to the initial ground state. For this asymptote coordinate choice, Q a v(q) is the overall av pseudofermion 
(or hole) phase shift whose value is given only to within addition of an arbitrary multiple of 2tt and 2ir <& a v. a *v' (<7j , qy ) 
is an elementary two-pseudofermion phase shift. 

The studies of Ref. [33j reveal that the overall pseudofermion phase-shift choice Q av (q) = Q° av + Q% v {q) associated 
with the asymptote condition (|39[) corresponds to a generalization of the conventional phase-shift definition previously 
used in the BA literature for the particular case of the n — 1 and m — initial ground state. (All the discussions and 
analysis presented below in this paper also apply to the phase-shift definition Q a u(q) = <3a„ + Q*j,(<7)j provided that 
the av phase shifts 8 av (q) = Q av (q)/2 are multiplied by two.) 



B. THE TWO-PSEUDOFERMION PHASE SHIFTS: BARE-MOMENTUM DEPENDENCE AND 

LEVINSON'S THEOREM 



The scattering av pseudofermion or hole phase-shift Q^ v {q)/2 given in Eq. (0} results from the set of elementary 
two-pseudofermion scattering events associated with the transition from the "in" state to the "out" state (excited 
energy eigenstate). In contrast, the transition from the ground state to the "in" state has no scattering character and 
leads to the scatter-less phase-shift Q a au /2 = 0, T 7T /2 of Eq. I|17fl . 

The bare-momentum two-pseudofermion phase shifts a > v i (q, q') in units of it are related to the rapidity two- 
pseudofermion phase shifts $ ai/i a >„> (r, r') by Eq. @). The latter phase shifts are defined by the integral equations 
(Al)-(A13) of Ref. |(|. In Appendix B we provide a set of simplified equations which define the rapidity elementary 
two-pseudofermion phase shifts "Iw, Q v (^ r') in the limit m — * 0. Furthermore, in that Appendix we also provide 
closed form expressions valid in the specific limit m — ► and n — > 1. 

In figures 1-6 we plot the two-pseudofermion phase shifts which contribute most to these properties as a function 
of the bare-momenta q and q' for electronic density n = 0.59, spin density m — > 0, and U/t — > , U/t = 0.3, 4.9, 100. 
(The electronic density value n = 0.59 and the value U/t — 4.9 are those used in Ref. [24j for the description of the 
TCNQ photoemission dispersions observed in the quasi- ID organic compound TTF-TCNQ.) Analytical expressions 
valid for U/t — > as m — > of the same two-pseudofermion phase shifts are provided in Appendix C. 

Note that the function & c o,ci(q, q') plotted in Fig. 5 can have values in the domain <E> c o iC i(g, q') 6 (—1, 1). Thus, 
within the standard quantum non-relativistic scattering theory phase shift definition given in Eq. (|36|) . the corre- 
sponding phase shift it $ c o, c i(q, q') has values in the domain it <fr c o tC i(q,q') € (— it, tt). Note that the width of this 
domain is 2tt, whereas the definition (|39|l would lead to a domain width of 4ir. 

For m — > the two-pseudofermion phase shifts have a similar qualitative behavior for a large domain of U/t values 
corresponding to U/t > 4. Thus, in order to obtain further information about the q and q' dependence of such phase 
shifts, it is useful to derive an analytical expansion mt/U for m — > 0, which provides the exact behavior for U/t^$> 1. 
By use of manipulations similar to those reported in Appendix C, we find that for m — ► and U/t 3> 1 the elementary 
two-pseudofermion phase shifts plotted in Figs. 1-6 can be written as follows, 



7T$ 



cO, cO 



(9, q') 



it (Co ~ l/£o) / sin(g) - sin(g / ) 
2 V 2 sin(Trn) 



(Co + l/Co) ,1 sin(g') 



sin(7m) 



(40) 



t 3>co si(q, q') = — arc tan ( sinh ( sin(q) + arc sinh ( tan ( — 

2 V V U V V n 



(Co - 1) 2q< 



(41) 



7r<i>si,co(<7, q) = —- arc tan ^sinh ^arcsinh ^tan^ — 
sgn(q)7r 



2a/2 



q = ±k F , 



-jj- sm (q ) 



q / ±k F 



(42) 
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FIG. 1: The elementary two-pseudofermion phase shift "I> c o, c o(q, q) in units of n as a function of q and q' for n = 0.59, m = 0, 
and (a) U/t — » 0, (b) U/t = 0.3, (c) C7/ 1 = 4.9, and (d) U/t = 100. The bare-momentum values q and q correspond to the 
right and left axis of the figures, respectively. 
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FIG. 2: The elementary two-pseudofermion phase shift "I> c o, s i(q, q) in units of 7r as a function of q and q' for n = 0.59, m = 0, 
and (a) (7/i — » 0, (b) U/t = 0.3, (c) (7/t = 4.9, and (d) U/t = 100. The bare-momentum values q and q correspond to the 
right and left axis of the figures, respectively. 
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FIG. 3: The elementary two-pseudofermion phase shift & 3 i t c o(q, q) in units of 7r as a function of q and q' for n = 0.59, m = 0, 
and (a) U/t — » 0, (b) U/t = 0.3, (c) (7/t = 4.9, and (d) U/t = 100. The bare-momentum values q and q correspond to the 
right and left axis of the figures, respectively. 
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FIG. 4: The elementary two-pseudofermion phase shift & s i t 3 i(q, q) in units of 7r as a function of q and q' for n = 0.59, m = 0, 
and (a) U/t — » 0, (b) = 0.3, (c) U/t = 4.9, and (d) U/t = 100. The bare-momentum values q and q' correspond to the 
right and left axis of the figures, respectively. 
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FIG. 5: The elementary two-pseudofermion phase shift $ c o, ci(q, q') in units of n as a function of q and q' for n = 0.59, m = 0, 
and (a) (7/i — » 0, (b) U/t = 0.3, (c) (7/t = 4.9, and (d) U/t = 100. The bare-momentum values q and q correspond to the 
right and left axis of the figures, respectively. 
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FIG. 6: The elementary two-pseudofermion phase shift & 3 i t c i(q, q) in units of 7r as a function of q and q' for n = 0.59, m = 0, 
and (a) (7/i — » 0, (b) U/t = 0.3, (c) (7/t = 4.9, and (d) U/t = 100. The bare-momentum values q and q correspond to the 
right and left axis of the figures, respectively. 
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P OC 

JO 



oo sin I a; ^ [arc sinh I tan I £ j ) — arc sinh I tan I ^ 



i sin(Trn) 2q' fq\ 

77 k cosl- ; q^±k F 

U 2 n \nJ 



sgn(g)7r 
2^2 

= [sgn(g)]7r(-i= - l) ; g = g' = ±fc F , , (43) 



g = ±kp , g' 7^ ztkp 
3 

-2^1 



7r$ c0 , c i( 9 , g') =arctan(-^sin(< ? )+tan(|^) + - X) | , (44) 

and 

7T$ s i c i(g, g') = 4 sin(7T(5)^ cos(-) , (45) 

respectively. In these equations d = 1 — n and £0 is the interaction-dependent parameter which for zero spin density 
appears in the expressions of the low-energy quantities 0, IE El El ES IS EH EM EE E3- For instance, it is 
defined in Eq. (74) of Ref. [33. In the above U/t ^> 1 two-pscudofcrmion expansions ijlOjl. J31J, and (l4~4l . £ should 
be replaced by its first-order expansion in t/U, 

4t 

£0 = 1 + -T7 In 2 sin(Trn) ; U/t > 1 . (46) 

Note that the t/J7 leading-order term of the quantity (£0 + 1/Co) appearing in Eq. (|40|) is of second order. However, 
since the t/U second order term of the parameter £0, which is not given in Eq. 146(1 . does not contribute to (£0 + l/£o), 
that leading-order term should be considered. 

For U/t = 100 and n = 0.59 the bare-momentum dependence of the two-pseudofermion phase-shift analytical 
expansions (|40|I - H45|I is similar to the exact dependence plotted in Figs. 1-6. In figure 7 such phase-shift expansions 
are plotted for U/t — 4.9 and n — 0.59. In spite that they reproduce the exact behavior for U/t » 1 only, note that 
their use for U/t = 4.9 is indeed a reasonably good approximation. This can be confirmed by comparison of the two- 
pseudofermion phase-shift expansions plotted in Fig. 7 with the corresponding exact phase-shift values plotted in Figs. 
1-6. This result is consistent with the U/t >> 1 physics being dominant from U/t « 4 until U/t—* 00. It follows that 
the phase-shift analytical expressions (|4L)|> - (|45|) can be used as an approximation for the two-pseudofermion phase-shift 
expressions for such a range of U/t values, becoming an increasingly better approximation as U/t increases. 

In order to confirm that our theory is consistent with the expected general properties of the standard scattering 
theory, let us check whether the phase shifts it $ a!/j a'v'{q, q') associated with the elementary two-pseudofermion 
scattering events obey Levinson's Theorem |40j |. Such a theorem just states that when in the reference frame of the 
scattering center the momentum of the scatterer tends to zero the phase shift is given by 7riV& , where Nb is the number 
of bound states. In that frame the phase shift it <& au . a 'v'(q, q') reads it $ QI/j a 'u'{q — g', 0). Moreover, in our case there 
are no bound states and thus Nb = 0. Therefore, for the pseudofermion scattering theory Levinson's Theorem requires 
that, 

lim 7r$ a „, aV (g-g',0) = 0. (47) 

q-q'->0 

The validity of this result is confirmed by suitable analysis of the integral equations (Al)-(A13) of Ref. @, which 
reveals that tt $ Ql/: a i v i (r, r') = —n$> a v,a'i>' (— r,— r'). This result combined with the use of Eq. © and the odd 
character of the ground-state rapidity functions, such that A^ v (q) — — A"„(— q) [l5|, leads then to, 

TT <&cw, a'v' (q, q) = ®au, a'v' (— <Z, -?') ■ (48) 

This latter symmetry implies that 7r <& a i>, a'v'(q — g',0) is a odd function of g — g', what confirms the validity of the 
Levinson's Theorem (|47|) . 

Finally, let us address the relation of the pseudofermion phase shifts to the phase shifts considered in Refs. |38l l39]| . 
It is straightforward to show that the elementary phase shifts defined by Eqs. (32)-(38) of Ref. [3^ correspond to a 
particular case of the rapidity elementary two-pseudofermion phase shifts defined by the integral equations (A1)-(A13) 
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FIG. 7: The U/t » 1 two-pseudofermion phase-shift expansions I|4()|l - I45|| in units of n for (a) $ c o, co(q, q), (b) $ c o, si(g, q'), 
(c) $ B i, c o(g, g'), (d) <E> s i, 3 i(g, g'), (e) $ c o,ci(g, g')> and (f) 4>si, c o(g, q), respectively, as a function of q and q' for n = 0.59 and 
{//t = 4.9. The bare-momentum values q and q' correspond to the right and left axis of the figures, respectively. For U/t — 4.9 
the phase-shift large U/t expansions H4UII-I45II are a reasonably good approximation for the corresponding exact phase shifts 
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of Ref. 0. If one considers the PS subspace spanned by energy eigenstates such that N av = for the av ^ cO, si 
branches and L a _i/ 2 = for a = c, s, the general integral equations (A1)-(A13) of Ref. M reduce to the integral 
equations (32)- (38) of Ref. Thus, the phase shifts previously considered in Refs. [38l l39l | correspond to a 

particular case of the general elementary two-pseudofermion phase shifts. In contrast to the interpretation of Refs. 
|38ll3^ |. the scatterers and scattering centers associated with the phase shifts defined by Eqs. (32)-(38) of Ref. [3^ are 
the pseudofermions and pseudofermion holes, rather than the corresponding pseudoparticles and pseudoparticle holes 
considered in that reference. Indeed, under the ground-state - excited-state transitions the discrete bare-momentum 
values of the pseudoparticles do not acquire the scattering momentum shift Q% u {qj)/L given in Eq. QJ. Instead, 
within the pseudoparticle representation the functional Q% u {q_j) appears in the energy spectrum, where it leads to 
/-function energy terms associated with the two-pseudoparticle residual interactions 6]. In turn, such interactions 
do not exist for the pseudofermion representation. 



C. INVARIANCE UNDER THE ELECTRON - ROTATED-ELECTRON UNITARY TRANSFORMATION 

In the previous section it was mentioned that the pseudofermions are not in general invariant under the electron 
- rotated-electron unitary transformation. However, the exception is for the av ^ cO, si branches as the canonical 
momentum q approaches the limiting values, q —> ±g°„. This is consistent with the result of Ref. 5] that creation of 
one cv ^ cO pseudoparticle at iq^, = ±[ir — 2fci?] (and one sv ^ si pseudoparticle at ±q°„ = ±[fc_F \ — fcp jj) leads to a 
change v in the number of lattice sites doubly occupied by both electrons and rotated electrons (and singly occupied 
by both spin-down electrons and spin-down rotated electrons). The same result applies to the corresponding cv ^ cO 
pseudofermion (and sv ^ si pseudofermion). Here we study the role of such a symmetry in the scattering properties 
of the av ^ cO, si pseudofermions for initial ground states with densities in the ranges < n < 1 and < m < n. 
(As mentioned above, the problem is addressed for initial ground states with densities n = 1 and/or m = in 
Sec. IV-D.) First we use the phase-shift expressions involving av ^ cO, si pseudofermions of canonical momentum 
q = ±q°„ to find important information about their properties as scatterers and scattering centers. Next we use 
such properties to clarify the effects of their invariance under the electron - rotated-electron unitary transformation. 
Below in this subsection we use in general the notations av and a'v' to denote the branch of the pseudofermion or 
pseudofermion-hole scatterers and scattering centers, respectively. 

Our first goal is to show that the av ^ cO, si pseudofermions of limiting canonical momentum q — ±q av are 
not active scatterers. By active scatterers we mean scatterers whose overall phase shifts generated by the ground- 
state - excited-energy-eigenstate transitions lead to a shift of the corresponding canonical-momentum values. Let 
us confirm that the ground-state limiting canonical-momentum values q — of the pseudofermions belonging to 

av =^ cO, si branches are not shifted by the ground-state - excited-energy-eigenstate transitions. In contrast to the 
cO pseudoparticles and the usual band particles and Fermi-liquid quasi-particles, the band bare-momentum limiting 
values ±q°„ = ±7r [N* v — 1]/L of the av ^ cO, si pseudoparticles can be changed by the above transitions. Since 
±Aq°„ = ±7r AN* V /L, such an exotic behavior occurs when the deviation AN* U , Eq. <|18(l . generated by the ground- 
state - excited-energy-eigenstate transition is finite. From use of Eqs. (A8)-(A13) of Ref. Q we find that the two- 
pseudofermion shifts 7r $0,^ a i v i (l q®, u , q) which can be written as ir <& av ^ c „< (l[tt — 2kp ], q) and ir $ QI/ , sv i {t[kp^ — kp jj, q) 
for av = cv =/= cO and av = sv si, respectively, have the following expression, 

n$ a v,<x'v>{Lq av ,q) = -y <W',co(£<*,c - S a , s ) + +v + v' - \v - v'\) ; i = ±1 ; av ^ cO, si , (49) 

where q ^ i-Qav f° r a ' v ' = au ' 7^ si. Use of this two-pseudofermion phase-shift expression in the general overall 
scattering phase-shift expression (0} leads to, 
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AN au -[S a , c -5 a , s ]AN c0 -y"(v + v' -\v-v'\)AN av , ; l = ±1 ; av ^ cO, si . (50) 



Comparison of this expression with that of iAq al , = lit AN*^/ L where N*^ — N aL , + N^ v and N^ v is provided 
in Eq. l|2~4T confirms that Q% u {^ v )/L = bAq au . The value of the two-pseudofermion phase shift P^l is not well 
defined for q ^ i^q^v and a'v' — av' ^ cO, si. However, from the rapidity expression ^aAq) = A^ tl/ (q(q)) for the 
PS excited energy-egenstates, where N° av (q) is the expression for the initial ground state [y], one confirms that the 
relation Q^^^q^) / L = iAq av is valid for all PS excited states. Indeed, that the rapidity functions of the excited 
energy states of a given initial ground state equal those of the latter state provided that in the argument of such 
functions the ground-state bare momentum is replaced by the excited-state canonical momentum implies that the 
corresponding bare-momentum and canonical-momentum bands have precisely the same width. 
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This result has a deep physical meaning: the scattering phase shift leads to a canonical-momentum shift 
Q^ v (±q au )/L that exactly cancels the bare-momentum shift ±Aq al/ . This implies that the overall canonical- 
momentum shift ±Aq av = ±Aq al/ + Q^ v (±q av ) / L indeed vanishes. It follows that for the av ^ cO, si pseudofermions 
the limiting canonical momenta have the same values, ±?„ w , both for the ground state and excited energy eigenstates 
and thus for q = ±q au such objects are not active scatterers. 

Let us next investigate the properties of such pseudofermions as scattering centers. By use of Eqs. (Al)-(A13) 
of Ref. @ and Eqs. (A.11)-(A.14) of Ref. O we find after some algebra that the two-pseudofermion shifts 
7r 4> ai/j a i v i (q, l q a i v i) which can be written as ir <I> Q1/ c „/ (q, l[tt — 2kp}) and ir <& at/j sv t (q, i[kp^ — kpi)) for a'v' = cv' ^ cO 
sv' ^ si, respectively, have the following general form, 



■K$av,aiv>{q,Lqa> v >) = ^ [S au> C Q - 8 a , a > (f + v' — \v - v' \ ) 



+ ~2 ^2 <''\n$avco(q,t''2kF) - 5 a ',s2n$ al , s i(q,i' kpi) ; t = ±l; a'v' ^ cO, si . (51) 
t '=±i 

Here av is any of the branches with finite pseudofermion occupancy for the excited state under consideration and the 
values of q are such that \q\ < q a , v , for av — a'v ^ cO, si and otherwise can have any value and thus correspond to 
all active av scatterers of that state. The form of these exact two-pseudofermion phase-shift expressions reveals that, 
except for the constant phase-shift terms, creation of one cv' ^ cO pseudofermion (and one sv' ^ si pseudofermion) at 
canonical momentum l[it — 2k f] (and L[kp-\ — kpi]) is felt by the av pseudofermion or hole active scatterers as a shift 
m/L of both cO bare- momentum Fermi points (and a shift m/L of both cO bare- momentum Fermi points and a shift 
—l2tt/L of both si bare-momentum Fermi points). Thus, such scatterers effectively feel that they are scattered by 
cO fermi-point current shifts (and cO and si fermi-point current shifts), rather than by the cv' ^ cO (and sv' ^ si) 
pseudofermion created at canonical momentum l[tt — 2kp] (and i[kp-\ — ^fj,])- 

Active scattering centers are those which contribute to the scattering phase shift (0J. For instance, small- momentum 
and low-energy cO and si pseudofermion particle-hole processes in the vicinity of the Fermi points, called elementary 
processes (C) in Ref. 0, do not generate active scattering centers. Indeed, within such processes the phase shifts 
generated by the pseudofermion particle excitations exactly cancel those produced by the pseudofermion hole excita- 
tions. The part of the bare-momentum distribution-function deviation generated by a'v' active scattering centers can 
be written as AN^, F , (q') + AN F U , (q'). Here AN^, F ,(q') is generated by the processes called elementary processes (A) 
in Ref. 0], which create and annihilate (and create) a'v' — cO, si pseudofermions (and a'v' ^ cO, si pseudofermions) 
away from the Fermi points (and away from the limiting values ±g° v ,). In turn, AN F , v ,(q') is generated by the pro- 
cesses called elementary processes (B) in the same reference, which create and annihilate (and create) a'v' — cO, si 
pseudofermions (and a'v' ^ cO, si pseudofermions) at the Fermi points (and at the limiting values ztq^,^, ). In this 
subsection we are mostly interested in the scattering centers associated with the deviation AN F , u ,(q'), whose general 
expression reads, 



AN F 



n f , 



3q',l,q° , , i a ' V ' = fi l 



<V, <.«<>, , ; QtV^cO, si. (52) 



Here the deviation numbers (and numbers) AN F , V , (and N F , ,) are such that AN a /^ = AN F , V , + AN^, F (and 
AW = N F V ,+N% F ). They can be expressed as AN F U , = AN F v , t +1 +AN F U ,^ (and N F V , = A^.+i+^'.-i), 
where AN F , V , ±1 is the deviation in the number of a'v' pseudofermions at the right (+1) and left right (—1) Fermi 
point (and N F , U , is the number of a'v' pseudofermions created at iq a i v i with l = ±1). The associated deviation 
current numbers (and current numbers) read 2AJ F , V , — AN F , V , +1 — AJV^, _ 1 (and 2J F , U , = N F , V , +1 — N F , V , 

According to the PDT of Rcfs. the elementary processes (A), (B), and (C) mentioned above lead to qual- 

itatively different contributions to the spectral- weight distributions. The PDT studies of these references considered 
that creation of cv' ^ cO and sv' ^ si pseudofermions at the limiting bare-momentum values is felt by both the cO 
and si scatterers as effective cO scattering centers and effective cO and si scattering centers, respectively. However, 
in these studies that was only considered for cO and si scatterers of bare-momentum value q = ztq Fa ^. We emphasize 
that the general two-pseudofermion expression l|51|) generalizes that result to all active av scatterers, including both 
av = cO, si and av ^ cO, si active scatterers of arbitrary bare momentum q. 

£ From the linearity in the deviations of the overall scattering phase shift 10} one can write Q^ y {q)/2 = Q%i NF \q)/2+ 

Qa» F \q)/2 where Qai NF \q)/2 and Qai F \q)/2 result from the a'v' scattering centers associated with the deviations 
AN^, F (q') and AN F u ,(q'), respectively. Also the part of the total momentum deviation (|2*U|l associated with the 
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elementary processes (A) and (B) can be written as AP NF + AP F . After some algebra involving the use of Eqs. 
i|5T)l. and (|52"|> we reach the following expressions for these quantities, 
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(53) 



— 2 — = 2^ 2^ "^M'^?! 1 ?*w) — ^ — 

qV-cO, si U — ±1 



+ ^ t'TT $ QV , co(?, t' 2fc F ) [A Jf + J cu> + E J < 

(,'=±1 J/'=l 

oo 

+ ^ i' 7T $ a „, ,1 (q, L 1 k Fi ) \ A J F - 2 ^ J F SV , 



and 
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</=±l v'=1 

^ cO - #a, a' (v + l/ -\U- v'\) 

a' v'^cO, si 



E E?' A <^') + E EK-fl^^W); 

aV— cO, si/ 7 g' cf'^cO 



u a!v' i 



(54) 



AF F = tt[L 



c, -1/2 



53 + 4fc F [ a jf + 53 j f , + 53 + 2k Fl [ a j f - 2 53 J F , 



(55) 



i/'=i 



j/'=2 



v'=2 



The general expression of the phase shift Qai F \q)/2 given in Eq. (|54|l is valid for all active av scatterers which 
for the av ^ cO, si branches refer to bare-momentum values in the range \q\ < q® av . In the AP F expression of Eq. 
(155(1 we have included the contribution from the —1/2 Yang holons. (The momentum contributions from the +1/2 
Yang holons and ±1/2 HL spinons vanish.) Note that the current contributions to the momentum spectrum AP F 
given in Eq. I|55f). which multiply twice the value of the cO and si Fermi momenta 2kp and kp^ are identical to 
the current contributions to the scattering phase shift l(54|l which multiply the phase shifts 7r <fr alJ , c o(q, 2kp) and 
7T <& a v, si (<?, kpi ) , respectively. 

That the cv' ^ cO pseudofermions and sv' ^ si pseudofermions created at limiting canonical-momentum values are 
not active scatterers and are felt by the active scatterers as effective cO scattering centers and cO and si scattering 
centers, respectively, follows from their invariance under the electron - rotated-electron unitary transformation. Such 
properties reflect the following important decoupling: as q — > ±q^, v ,, the cv' ^ cO pseudofcrmion (and sv' ^ si 
pseudofermion) separates into 2v' independent holons (and 2v' independent spinons) and a cv' (and sv') FP current 
scattering center. By independent holons and spinons we mean those which remain invariant under the electron 
- rotated-electron unitary transformation. It follows that the Yang holons and HL spinons are also independent 
holons and spinons, respectively. By a cv' (and sv') FP current scattering center we mean the elementary current 
J F V , — l/2 (and jf v , = t/2) generated by creation of one cv' + 1 cO pseudofermion (and one sv' ^ si pseudofermion) 
at q = l[tt — 2k F ] (and q = i[kp^ — kpi\). As confirmed by the form of the phase shifts given in Eqs. (|5*H and J^U, 
such elementary currents are felt by the av active scatterers as elementary shifts of both cO Fermi-points (and both 
cO and both si Fermi-points). This justifies the designation FP (from Fermi-points) current scattering center. 

The cv' cO and sv' ^ si pseudofermions of limiting canonical-momentum value ±q^, v , have energy 2v' /i and 
2v' [1q H, respectively, relative to the ground state. The energy of one —1/2 independent holon and one +1/2 indepen- 
dent holon (and one —1/2 independent spinon and one +1/2 independent spinon) relative to that state is 2\x and zero 
(and 2/i H and zero), respectively. Thus, for q — ±q®>v> the cv' pseudofermion (and sv' pseudofermion) energy is 
additive in those of the corresponding v' —1/2 holons and v' +1/2 holons (and v' —1/2 spinons and v' +1/2 spinons). 
This is consistent with the a'v' ^ cO, si pseudofermions loosing their composite character as q — > ±q° V /. While the 
whole pseudofermion energy goes over to the v' —1/2 holons or v' —1/2 spinons, part or the whole pseudofermion 
momentum, respectively, is transferred over to the FP current scattering center. In contrast, for \q\ < q^i u , the 
cv' pseudofermion (and sv' pseudofermion) energy relative to the ground state reads e cl / (q) = 2v' p, + e° cv , (q) (and 
e sv .{q) = 2v'hqH + e° sv ,(q)), where the energy e° a , u ,(q) ^ is defined in Eqs. (C.17), (C.18), (C.20), and (C.21) 
of Ref. Q. That the energy e° a i v i(q) is finite is consistent with the composite character of the 2v' holons and 2v' 
spinons of the cv' and sv' pseudofermion, respectively. In this general case the pseudofermions are not invariant under 
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the electron - rotated-electron unitary transformation. The holon or spinon degrees of freedom are then combined 
with the pseudofermion scattering part and, therefore, the cv' 7^ cO and sv' 7^ si pseudofermions created under a 
ground-state - excited-state transition are felt by the av pseudofermion and av pseudofermion-hole active scatterers 
as independent scattering centers, unrelated to the cO and si Fermi points. 

In what the independent holons associated with a cv' pseudofermion of canonical momentum q — iq^, = ±[ir— 2kp] 
is concerned, the main difference between having (i) L c _ 1 / 2 = v' —1/2 Yang holons plus £ c ,+i/2 = v ' +1/2 Yang 
holons and (ii) one cv' pseudofermion is that the 2v' Yang holons have a total ?7-spin value S c — v' whereas for the 
2v' holons associated with the pseudofermion the total 77-spin value is S c = and thus corresponds to a 77-spin singlet 
configuration. (The same holds for the sv' pseudofermions provided that one replaces holons by spinons and Yang 
holons by HL spinons.) As the Yang holons and HL spinons, the 2v' independent holons (a' = c) or 2v' independent 
spinons {a' = s) of a a'v' pseudofermion of canonical momentum q = ztq 1 ^,^, are neither scatterers nor scattering 
centers. 

Within the PDT of Refs. 0, the av = cO, si canonical-momentum Fermi-point deviations tAqp av play 
a central role in the spectral-function expressions through the related quantity 2A' aL , = [t Aq~Fuv / {2tt / L)] 2 where 
iAqFav/ (2n/L) — lAN^ v t + Q"% v {iq Fav )l2-n . We thus emphasize that the validity of the corresponding expressions 

(37) and (40) of Ref. 0] follows from the general expression for the phase shift Q%i F \q)/2 given in Eq. 

D. INVARIANCE UNDER BOTH THE ELECTRON - ROTATED-ELECTRON AND 
PSEUDOPARTICLE - PSEUDOFERMION TRANSFORMATIONS FOR n = 1 AND m = DENSITIES 

The general pseudofermion scattering theory also applies to initial ground states with densities n = 1 and/or m — 0, 
provided that the specific features reported here are taken into account. For an initial ground state with electronic 
density n — 1 (and spin density m — 0) one has that N*„ = (and N* u = 0) for the cv (and sv) band and thus 
the corresponding pseudofermion branch does not exist. For simplicity, we focus our attention onto excited energy 
eigenstates of such an initial ground state with a single cv 7^ cO pseudofermion (and a single sv 7^ si pseudofermion). 
For these excited states, N* v = 1 (and iV*„ = 1) and the corresponding cv (and sv) bare-momentum band reduces to 
the bare momentum zero. 

A property specific to n = 1 (and m = 0) initial ground states is that a transition to an excited state involving 
creation of one cv 7^ cO pseudofermion (and one sv 7^ si pseudofermion) always also involves creation of 2v cO 
pseudofermion holes (and 2(v — 1) si pseudofermion holes). Furthermore, if one considers the n = 1 and m = initial 
ground state it follows from Eqs. (|B"§|) . I|B21|I . (|1522|) . and (|1523[) of Appendix B that $ c0 , sv > (q, 0) = $ C!/) sv > (q, 0) = 
®su. cv'iq, 0) = for v > 1 and v' > 1. In turn, we find below that the value of the phase shifts it $ c o, C v(q, 0) (and 
7T $si, si/(?, 0)) is fully determined by the 2v (and 2(y — 1)) bare-momentum values of the cO pseudofermion-hole (and 
2(y — 1) si pseudofermion-hole) scattering centers. Thus, for excited states of the n = 1 and m — initial ground 
state the cO and si scatterers feel the created cv ^ cO pseudofermion (and sv ^ si pseudofermion) as cO effective 
scattering centers (and si effective scattering centers). 

These effective scattering centers are different from those considered above for the excited states of ground states 
with electronic density (and spin density) in the range < n < 1 (and < m < n). Indeed, for the excited states of 
a n = 1 (and m = 0) initial ground state considered here the current number J^ v (and Jf v ) vanishes and thus there 
are no cv 7^ cO (and sv 7^ si) FP current scattering centers. The occurrence of the type of cO (and si) effective 
scattering centers considered in this section follows from the non-scatterer character of the corresponding cv 7^ cO 
pseudofermions (and sv 7^ si pseudofermion), as discussed below. 

In spite of the lack of ground-state cv 7^ cO (and sv 7^ si) pseudofermion bands, the scattering theory can be 
generalized to an initial ground state with electronic density n = 1 (and m — 0): we recall that the "in" asymptote 
one-pseudofcrmion scattering states do not contribute to the direct-product expression of the initial ground state but 
rather to that of the "in" state defined in Sec. III-C. For the excited states of n — 1 (and m = 0) ground states 
considered here, the cv 7^ cO (and sv 7^ si) bare-momentum band corresponds to a single value at q = 0. For initial 
ground states with densities < n < 1 and < m < n the scattering canonical-momentum shift Q* ly (g)/L, Eq. 
has the same value whether one uses the ground-state rapidity functions A° w (g) and hP a , v ,{q') or the corresponding 
"out" -state (excited-energy-eigenstate) rapidity functions h au (q) and A Q /„/(g') in the general expression |JB} for the 
two-pseudofermion phase shifts. Indeed, these two alternative definitions of the two-pscudofcrmion phase shifts lead 
to the same value for the functional Q^(q)/L up to contributions of order 1/L. In turn, the general pseudofermion 
scattering theory introduced above also applies to initial ground states with densities n = 1 and/or m = provided 
that the following procedure is performed: 

- Since for a n = 1 and/or m = initial ground state there are no cv 7^ cO and/or sv 7^ si pseudofermion bands, 
the two-pseudofermion expression |JBJ should be replaced by $ a!/) a i v i (q, q') = §> av ^ a i u i(4:tA au (q)/U,4tA a > v >(q')/U), 
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where the rapidity function A au (q) (and A a i v r(q')) is that of the excited state and av — cv with v > and/or av = sv 
with v > 1 (and aV = cv' with i/' > and/or a'u' = sv' with v' > 1), and otherwise is that of the initial ground 
state. Since the former rapidity functions are those of the excited state under consideration, it follows that for the 
particular case of such an initial ground state the quantity © is a functional rather than a function. 

We found in the previous subsection that for initial ground states with densities in the ranges < n < 1 and 
< m < n, the av ^ cO, si pseudofermions with limiting canonical momentum given by ±g^ y are not active 
scatterers. As discussed below, for the excited energy eigenstates considered in this subsection the cv ^ cO (and 
sv si) band reduces to a single discrete canonical momentum value at q = 0. Such excited states have electronic 
density n — > 1 (and spin density m — ► 0) and thus the single q = value corresponds to the limiting canonical 
momentum values ±gj?„ = ±[ir — 2kp] (and ±q a sv = ±[fcpf — fcpj) such that ±(f{?„ — > as n — > 1 (and ±q a sl/ — > 
as m 0). It follows that similarly to the pseudofermions considered in the previous subsection, the cv ^= cO 
pseudofcrmion (and sv ^ si pseudofermion) of the excited states considered here are invariant under the electron - 
rotated-electron unitary transformation. Indeed, it follows from the results of Ref. p| that creation of such a cv ^ cO 
pseudofermion (and sv ^ si pseudofcrmion) leads to a change v in the number of lattice sites doubly occupied 
by both electrons and rotated electrons (and singly occupied by both spin-down electrons and spin-down rotated 
electrons). As a result, they separate into 2v independent holons (and spinons). Such 2v independent holons (and 2v 
independent spinons) are fully decoupled from the above mentioned cO effective scattering centers (and si effective 
scattering centers) also associated with the cv ^ cO pseudofermion (and sv ^ si pseudofermion). It is then required 
that they are not active scatterers. However, for the class of excited states considered here a necessary condition for 
such objects not being active scatterers is that they are not scatterers at all. Indeed, once the cv ^ cO (and sv ^ si) 
bare-momentum band of the "in" 

state has a single value at q = 0, it is required that the corresponding canonical-momentum band of the "out" state 
has also a single value at q = 0. This implies both that <3*„(0)/Z/ = and Q a ^(0)/L = 0, and thus that the cv ^ cO 
(and sv =/= si) pseudoparticle remains invariant under the pseudoparticle - pseudofermion unitary transformation. 
Indeed, since the cv ^ cO (and sv =/= si) band does not exist for the initial ground state, one has that the scatter- 
less overall phase shift vanishes, Q a av /2 — 0, and thus for the above reasoning the overall phase shift is such that 
Qav(0)/2 = Q* 1 y(0)/2 = 0. It follows that the corresponding cv ^ cO (and sv ^ si) pseudofermion is not a scatterer. 

However, the requirement that such a cv ^ cO (and sv ^ si) pseudofcrmion is not a scatterer and thus that 
Q%(0)/L = (and Q%(0)/L = 0) imposes a specific form to the corresponding two-peudofermion phase shifts 
^cv, a >u' (0, q') (and a > v > (0, q')). Since such objects are neither scatterers nor scattering centers, the quanti- 
ties 7T $ cl/j Q v (0, q') (and n $ ;5 ^ ! a ' u ' (0, q')) are not real two-peudofermion phase shifts: they are just effective two- 
peudofermion phase shifts whose values are such that the overall scattering phase shift Qf„(0)/2 (and Q%(0)/2) 
vanishes. Thus, there is no requirement that they do obey Levinson's Theorem. Also the phase shifts ir cu(q' , 0) 
and 7r § a 'v', sv (<?', 0) are not required to obey that theorem for aV = cO, si, once they refer to effective cO and si 
scattering centers, as confirmed below. 

Let us consider three types of excited energy eigenstates of the above class. Those are states with finite pseud- 
ofermion occupancy for the cO and si bands plus (a) one cv ^ cO pseudofermion and one sv' ^ si pseudofermion, (b) 
one cv ^ cO pseudofermion, and (c) one sv' ^ si pseudofermion. For simplicity, we consider excited energy eigen- 
states of the n = 1 and m = initial ground state. Such a ground state is described by full cO and si pseudofermion 
bands whose Fermi bare momentum reads 2kp = n and kp ^ — kp^ = kp — n/2, respectively. Thus, from the use 
of Eq. (B.ll) of Ref. H| we have that AA^ c0 = -AA^ = ~2v, AN s i = -v - v', AN^ = 2[y' - 1) for the excited 
states (a), AiV c0 = -AA^ l = -2v, AN sl = -v, and AN^ = for the excited states (b), and AiV c0 = -AA ^ = , 
AiV s i = —i/, and AN^ = 2(v' — 1) for the excited states (c). On the other hand, according to Eqs. (|B9|I . IjBlOp . 
and l|B23|) of Appendix B, for m — > and n — ► 1 the two-pseudofcrmion phase shifts that contribute to Q%(0)/2 

and Q*,(0)/2 simplify to # ci/ , 8V , (r, r') = $ su >,co(r',r) = ^ SV ',cu{r',r) = 0, § cw , o0 (r, r') = iarctan(^), and 

$s>v,si (r',r) — i arctan f p^l ) ^ or v ' > ^ an d v > 0. It follows that for the excited energy eigenstates (a)-(c) the 
equation Q% (0)/2 = and/or Q% ,(0)/2 = leads to, 



arctan 

2 = 1 

2(f'-l) 



[A cv (0) - A° c0 ( qi ) 



0, v > 0, 



arctan ( jy^)U i^' ( ° } " ) = ° ' > 1 ' ^'^ 



Here the first and second equations refer to the cO branch and both the states (a) and (b) and to the si branch 
and both the states (a) and (c), respectively. In these equations the set of 2v values {qi} and of 2(v' — 1) values 
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{q'l} correspond to the excited-energy-eigenstate cO pscudofermion holes and si pseudofermion holes, respectively. 
Moreover, the ground-state rapidity functions A® (q) = sinfco(g) and A° 1 (g) can be defined in terms of their inverse 
functions given in Eq. (A.l) of Ref. [15J. For n = 1 and m = the expressions provided in the latter equation lead 
to, 

, , s nfj sin(w sinfc (9)) T/ s [°° , sinfo; A° sl (q)) . . 

q=k (q)+2 du 2 j, 4 »l 9=/ 4^777^ J ° W ' 57 

Jo lu (1 + e- !ujU '^ t ) J oj cosh(oj (7 / 4i) 

where Jo(w) is a Bessel function. These equations define the inverse of the functions fco(<?) and A^g), respectively. 
(To arrive to the second equality of Eq. we used that A° 1 (±/cf) = A° 1 (±7r/2) = ±oo.) The form of the equalities 
given in Eq. (|Sl)jl reveals that the corresponding solutions A c „(0) = A cv (0, {qi}) and/or A siy i(0) — A sv > (0, {q\}) are 
functions of the above sets of bare-momentum values {qi} and {q' ; }, respectively. 

We emphasize that the solution of the BA equations (13)-(16) of Ref. Q for the above excited states leads to 
functions A CI/ (0) = A c „(0, {qi}) and A su /(0) = A su / (0, {q'i}) for the rapidities A CI/ (0) and A SI ,<(0), respectively, which 
also obey Eq. 156fl . This confirms that for such excited states the exact solution of the BA equations is indeed 
equivalent to imposing the symmetry requirement Q*„(0)/2 = and Qf u ,(0)/2 = associated with the non-scatterer 
character of the corresponding cv and sv' pseudofermion, respectively. 

The above functions A cly (0) = A cl/ (0, {qi}) and A s „/(0) = A si ,<(0, {q\}) are to be used in the following expressions, 

.*.,,. l( o,rt - , »„, ( 4t a a» , iii|M) _ , rctm ( « ^ a. v ■» - wi ) , (58) 

so that Q%(0)/2 = and/or Qf v ,(0)/2 = 0. The simplest case corresponds to v — 1 and/or z/ = 2 such that solution 
of Eq. @2i) leads to A cl (0,gi, 92 ) = [A° ( 9l ) + A°,(g 2 )]/2 and/or A s2 (0, 9 ' 2 ) = [A^^'O + A° 1 (g' 2 )]/2, respectively. 

The requirement that the cv ^ cO pseudofermion (and sv ^ si pseudofermion) considered here is not a scatterer 
implies that the corresponding rapidity function A c „(0) = A cl/ (0, {<#}) (and A s!/ /(0) = A s „< (0, {<?';})) does not in 
general vanishes and is the unique solution of the first (and second) equality of Eq. i|56[) . Combination of this result 
with the two-pseudofermion phase shift of Eqs. (|B21|) and 1|B22(1 of Appendix B reveals that the cO scatterer two- 
pseudofermion phase shift 7r ^>da, cu{l, 0) and the si scatterer two-pseudofermion phase shift 7r $ s i. s „/ (q, 0) are such 
that, 



7r$cO,ct/(?jUJ = 7r$ c0 , c ^^ — , — j = -aictan^ — — j 

7r« s i, SI /'(g,0] = 7T$ s i )S1/ /^ — , — ) = arctan^ (^ -\)U 



= ±-j=; q = ±k F . (59) 

We note that in addition to the cO or si scatterer bare-momentum q, the two-pseudofermion phase shifts provided 
in Eq. I|59|l are functions of the set of 2v bare-momentum values {qi} or 2(v' — 1) bare-momentum values {q\} of 
the 2v cO pseudofermion-hole scattering centers or 2(y' — 1) si pseudofermion-hole scattering centers, respectively, 
also created under the ground-state - excited-energy-eigenstate transition. As a result of the creation of the cv ^ cO 
(and sv' ^ si) pseudofermion, the cO (and si) scatterers acquire the phase shift it <& c o,cv(q,0) (and ir $ s i. sv i (q, 0)) 
whose value is fully controlled by the 2v (and 2(v' — 1)) bare-momentum values of the 2v (and 2(v' — 1)) cO (and si) 
pseudofermion-hole scattering centers. Thus, through the {qi} (and {q'i}) momentum dependence of the phase shift 
tt^cO, cv(q,0) (and ir <fr s i ySV i(q,0)), the cO (and si) scatterers feel the created cv ^ cO (and sv' ^ si) pseudofermion 
as 2v cO effective scattering centers (and 2(v' — 1) si effective scattering centers). 

Similar results are obtained for excited states of initial ground states of density n = 1 and/or m — with finite 
occupancy for a larger finite number of av pseudofermions belonging to several av ^ cO, si branches, except that 
the number of equations defining the the set of rapidities {Aq,^} is in general larger than above and each of these 
equations is more involved than the two equations given in Eq. (|56[) . Importantly, such av pseudofermions are also 
invariant under the electron - rotated-electron unitary transformation. 



V. CONCLUDING REMARKS 



In this paper we have shown that the non-perturbative and strongly correlated scattering problem described in terms 
of electrons by the Hamiltonian considerably simplifies in terms of the pseudofermions ofRefs. 0, US El EE El: 
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We have confirmed here that in terms of pseudofermion scattering the spectral and dynamical properties are controlled 
at all energy scales and for all values of the on-site electronic repulsion by two-pseudofermion zero- momentum forward 
scattering only. This agrees with the preliminary analysis of the problem of Ref. j^jj . The matrix elements between 
the ground state and excited energy eigenstates of the corresponding spectral functions can be expressed in terms 
of the pseudofermion anticommutators |l4 [TBI IT^| . Our results show that such anticommutators are controlled by 
two-pseudofermion zero-momentum forward scattering through the associated pseudofermion and hole S matrices. 

Our results have also clarified the relation between the pseudofermion scattering properties and symmetry. Specifi- 
cally, it was found for the metallic pase at finite spin density that the invariance under the electron - rotated-electron 
unitary transformation of the av ^ cO, si pseudofermions with limiting canonical momentum q = ±g°„ implies that 
such composite objects separate into 2v independent holons (a = c) or 2v independent spinons (a = s) and a current 
excitation. The latter excitation is felt by all av active scatterers of arbitrary momentum as a shift of both cO Fermi 
points and shifts of both cO and si Fermi points, respectively. The effects of the invariance under both the electron - 
rotated-electron unitary transformation and pseudoparticle - pseudofermion unitary transformation were also studied 
for the Mott-Hubbard insulator at zero spin density. The property that only the holons and spinons that remain 
invariant under such a transformation are allowed to exist as independent quantum objects which are not part of 
2^-holon and 2i/-spinon composite pseudofermions, respectively, is general and also applies to the Yang holons and 
HL spinons. Such objects are neither scatterers nor scattering centers. 

That our choice of scatterers and scattering centers profits from the transformation laws under the electron - rotated- 
electron unitary transformation justifies that all "in" and "out" states of the theory are excited energy eigenstates. All 
these states can be written as a direct product of "in" asymptote and "out" asymptote one-pseudofermion scattering 
states, respectively. Such a property combined with the simple form obtained for the pseudofermion and hole S 
matrices is behind the suitability of the pseudofermion representation for the study of the finite-energy spectral and 
dynamical properties 0, ITHl Hfij . The studies of Ref. |33l | clarify the relation of the pseudofermion phase shifts 
and S matrices introduced in this paper to the corresponding quantities of the conventional spinon-holon scattering 
theory of Refs. [3(], El E3- The choice of scatterers and scattering centers of the latter theory is also based 
on associating the Bethe-ansatz quantum numbers with quantum objects, yet it does not profit from the invariances 
under the electron - rotated-electron unitary transformation. The clarification in this paper of the relation between the 
pseudofermion scattering mechanisms and symmetry provides further useful information about the PDT microscopic 
processes [l5lE^|. 

Since the transport and dynamical properties Edl and other properties predicted by the ID Hubbard model were 
observed in low-dimensional complex materials and the investigations presented in Refs. pijL l24l |25| confirm 
that the PDT describes successfully the unusual finite-energy spectral features observed by angle-resolved photo- 
electron spectroscopy in quasi-lD organic metals, our results also contribute to the further understanding of the 
non-perturbative scattering mechanisms behind these properties. While the studies of this paper considered the ID 
Hubbard model, our results are of general nature for many integrable interacting problems |43j and therefore have 
wide applicability. 
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APPENDIX A: PSEUDOFERMION SCATTERING PROCESSES ENERGY CONSERVATION 

Both the scattering phase shift Q% v {qj)/2, Eq. JI), and the overall phase shift Q a u{<ii)/^ = Q Q av l^ + Q% v {lj)/'^i 
Eq. (|16|) , conserve the total energy. 

Let us start by confirming that for L >> 1 the above overall phase shift conserves the total energy. The cO and cl 
pseudofermion momentum distribution function deviations of Eq. can be written as, 

Mf m %) = AAfWfa) + AA/-Wfe) ; av = cO, si . (Al) 

Here AAf a „(q~j) = AN au (qj) and the deviations AAfi]) (q~j) and ANaJ (q~j) are associated with the ground-state - 
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virtual-state transition and the virtual-state - "out"-state transition, respectively. The latter deviation reads 
AA/$ (?) = KM + Qau(q)/L) - Km 



= -j- sgn(g) Q QI /(sgn {q)q Fau ) 5{q Fau -\q\); a = cO, si . (A2) 



Here g = g is the virtual-state canonical momentum value, the Fermi momentum q Fau of the cO and si bands is given 
in Eq. JSJ, and Q av (q)/2 is the phase shift l|16|). Use of Eq. (|A2I) in the energy spectrum (0 leads to the following 
energy spectrum to first order in the canonical-momentum distribution- function deviations, 



AE 



(2) 

cO, si 



olv— cO,sl qj — — q° 

-is e ■ 

t — ±1 QI/— cO-Sl 



t — ±1 Q'i^ — cO . S 1 



(A3) 



In order to confirm that the factor — X) t =±i S Q iy= c o si 1 Qav{^ q Fa v) e av(q Fa v)i which multiplies 1/L in the last term 



on the right-hand side of Eq. (|A3|) . vanishes we have used the symmetry e QI/ (g) = e alJ (— q) and Eq. (|10|) such that 
t av (q Fav ) = for av = cO, si. (We recall that the occupancies of other ctv ^ cO, si pseudofcrmion branches vanish 
for the ground state.) Thus, we have just found that the energy contribution of first order in the canonical-momentum 
distribution-function deviations originated by the collective excitation (2) vanishes. 

Note that the energy (|A3ll decouples into two contributions, corresponding to the scatter-less phase shift Q°„/2 and 
scattering phase shift Q% y {tq Fav ) /2. This confirms that both these phase shifts conserve the energy independently. 

Since the collective excitation (2) involves all virtual-state cO and si pseudofermions, we used similar procedures for 
the evaluation of the energy contributions of order larger than one in the canonical-momentum distribution-function 
deviations. We find that all such contributions also vanish and decouple into independent and vanishing contributions 
corresponding to the scatter- less and scattering phase shifts. It follows that both the overall phase shift <5 Q1 /((?j)/2 
and the phase shifts Q°„/2 and <9*^(Qj)/2 conserve the total energy. 



APPENDIX B: ELEMENTARY TWO-PSEUDOFERMION PHASE SHIFT EXPRESSIONS FOR m -> 

AND n < 1 AND FOR m -» AND n -> 1 

We start by considering the limit m — > 0. The rapidity two-pseudofermion phase shifts ir <fr a u, a'v 1 (r, r') are defined 
by the integral equations (A1)-(A13) of Ref. [(|. By Fourier transforming these equations after considering that 
B = oo and thus r° = At B/U = oo for finite values of U /t, we arrive to the following equations valid for m — * and 
U/t finite, 

/+ro 
dr " A(r-r")n$ c0 , c0 (r", r') , (Bl) 
-ro 

7r$ ctM1 (r, r 1 ) = -iarctan(sinh(|(r-r'))) + j ° dr" A(r - r") n § c0 , s i(r" , r') , (B2) 



7r$ 8 i )C o(r, r') = -- arctan^sinh^(r - r , f , . 



1 f +r \ „ ^ c0 , e0 (r",r') 

' dr j 1 r- ; r f= ±oo 



-ro cosh/ ^(r — r")J 



sgn(r)7r 

= -^-; r = ±oo, (B3) 
7r$ sl , sl (r,r') = B(r-r') + 7 dr" ^^'^^ r '\ ■ r ± ±QO 



4 

sgn(r)7r 

-2V2 



ro cosh I 77 (r — r 



(f( 

; r = ±00 , r 7^ r 



= [sgn(r)]( — ^= - lW; r = r' = ±00, (B4) 



7r$ c0 ,oi/(r, r') = -arctan(- — -) + / dr" A(r - r") it $ c o, c „(r", r') ; ^>0, 

7r$ c0 , S! /(r, r') = 0; z/>l, 



* f m 1 /" +r ° w " *®co,cAr",n) 

7r$gi )t 3/(r, r ) = - / dr ? r ; > . 



4 



>-o cosh^(r — r")J 



7r*si,si/(?", r') = arc tan ^ — y) i r 7^ ±°° 5 i> > 1 , 



7T 

±— ■= ; r = ±oo ; v>\ 
v2 



7T $ SI /, co 0, r') = 7T $ s!/] 0I y (r, r') = ; 2/ > 1 

- ( V — T \ 

n$sv,si (r,r ) = arctan - - - J ; 1/ > 1, 



\ f T — T \ / 7* — T \ 

7T sw / (r, r ) = - 6„ „/ (r - r ) - arctan — ■ — - — - - arctan — ■ — - , v, v > 
2 \v + v' — 2/ \ v + v J 

In the above expressions the function Q ViV /(x) is defined in Eq. (B.5) of Ref. (4J, 



sin(W) i r (2+*4j r ( 1 *4 



o 



"(1 + e 2 ") 2 rfi-ijWl + if)' 



w tf dr 7T J l + e 2u) 

and 

4f sin Q 



ro 



U 



where Q is the parameter defined by Eq. (A. 5) of Ref. [la] and r(x) is the usual T function. 
Moreover, the two-pseudofermion phase shifts, 

7r$ cw , c0 (r,r ) = arctan(——) -- / dr — ; /r _ r „ ; ^>0, 



Z/> / 7T 



1 1 f +r ° 7r$ n / fr" rM 
7r^,^(ry) = ie„ jtf ,(r-r')-- / dr" Z^ c %t'L ; ^' >0 ' 

2 7T J_ ro + (-^) 2 ] 



7T 3>ci/ si/' (»", »" ) = / dr : 77— ; V > (J , 
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remain as in Ref. |6j. Furthermore, we note that the four expressions (|Bl|) - (|B4ll with cO — c, si = s, r = x, and 
ro = xo are equivalent to expressions (A9)-(A12) of Ref. j^. (For the phase shifts l|B3|) and (|B4I) this equality refers 
to values of r such that r ^ oo.) The rapidity phase-shift expressions given here for m — > correspond to some of the 
two-pseudofermion phase shifts plotted in units of n in Figs. 1-6. 

Let us now consider the limit m — ► and n — > 1. Note that the expressions i|B6|) and (|B8|) - (|B11|> are valid for 
m — > and all values of n such that n < 1. Thus, they also apply to the limit n — > 1. On the other hand, note 
that for all other rapidity two-pseudofermion phase-shifts the expressions for both m — > and n — > 1 are obtained by 
considering r = in the above integral equations. Such a procedure leads to closed form analytical expressions for 
all rapidity two-pseudofermion phase shifts. Thus, for m — ► and n — > 1 we find, 

7r$cO,co(r, /) = -B(r - r') ; 7r$ c0 , s i(r, r') = -- arctan(sinh(|(r - r'))) , (B18) 

T^ai,co(n^') = -^arctan(sinh(^(r-r'))) ; r ^ ±oo 
sgn(r)7r 

r = ±oo , (1j19) 



7r$ s i, sl (r, r') = B(r-r')\ r ^ ±oo 
sgn(r)7r 



2^2 

r = ±oo , r 7^ r 



2V2 ' 



[sgn(r)](-^=-l)7r; r = r' = ±oo , (B20) 



7r$ c0 , CI /(r, /) = -arctan(?— ^) , n$ a i, cv (r, r') = ir$> c0>sv >(r, r') = 0, i/>0,i/>1, (B21) 

- f T — t' \ 

K$si,su(r, r) = arc tan ^ — — -J; r^ioo; i^>l, 

= ±^=; r = ±oo; z/ > 1 , (B22) 
v2 

7r<IWo (r,/) = arctan(?— ^) ; tt $ c „, C!/ / (r, r') = i 6^, ^ (r - r') ; tt $ c!/) s „> (r, r') = , ^>0, (B23) 



and the phase-shift expressions (|T39|1 - (|T3 1 1|) are n independent and then have the same expressions as for n < 1. The 
two-pseudofermion phase shift expressions derived here are used in the derivation of the equalities of Eq. (|56|) . 



APPENDIX C: THE TWO-PSEUDOFERMION PHASE SHIFTS FOR U/t -> 

Here we derive expressions for the m — > two-pseudofermion phase shifts given in Appendix B that are plotted in 
Figs. 1-6 in units of tt, 7r$ c0 , c o(g, <?'), 7T$ o0 , s i(<2 q'), ir$ sl , c0 (q 7 q'), Tr$ sl , sl (q, q'), 7r$ c0 , c i(g, q'), and tt $ s i, c i(<Z, q'), 
for the limit U/t — *■ 0. Moreover, we also provide large-t//i expressions for the rapidity functions involved in the 
evaluation of the two-pseudofermion phase-shift expansions (I4UI1 - <|45|) . 

The evaluation of the two-pdseudofermion expressions provided here and of the expansions given in Eqs. 14011 - 145|) 
involves the use of Eq. © where the ground-state rapidity functions A^Lfg') are defined in terms of their inverse 
functions, whose expressions are given in Eqs. (A.l) and (A.2) of Ref. First, we use the latter expressions to 

derive the following closed- form expressions for the ground-state functions k°(q), A® (q), A®„(q), and A° 1 (g), valid for 
zero spin density m — > 0, values of the electronic density < n < 1, and limiting on-site repulsion values U/t — > 
and U/t > 1, 

k°(q) = |; \q\ <2k F , U/t - 

= sgn(g) [\q\ - k F ] ; 2k F < \q\ < ir/a , U/t -f 

= sgn(g)7r; \q\ = tt , U/i -> 

Atn 

= q - — ln(2) sin(g) ; |g| < tt ,[//*» 1 , (CI) 
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A° (g) = sin(|) ; |g| < 2k F , U/t^O 

= sgn(g) sin((|(z| - k F )) ] 2k F < \q\ < ir , U/t ^ 



0; 



\q\=n, 
2tn 



sin(g) - — ln(2) sin(2 9 ) ; |g| < tt , C//t » 1 , 



(C2) 



for v > 0, and 



A^(g) = sgnfo) sin( (|9| + 7rn) ) ; < |«| < (tt - 2fc F ) , C//t 

= ; g = 0, U/t^0 

= ±oo; g = ±(tt- 2fc F ) , U/t^0 

= ^ tan (^); o< | q | < (^-2fc F ), tvt»i, 



A° x (g) = sin(g); | 9 | <k F , U/t^0 
= ±oo; q = ±k F , U/t->Q 
U 



arcsinh) tan(— ) ) : \q\ < k F , U/t^>l, 

2nt \ n I 



(C3) 



(C4) 



respectively. 

Next, we use Eqs. I|C1|) - (|C4|) in the integral equations given in the Appendix A of Ref. |f|, which define the 
two-pseudofermion phase shifts. (Such equations were also presented in Ref. [44( with a slightly different notation.) 
By manipulation of these equations for the limit U/t — ► 0, we find the expressions for the above bare- momentum 
two-pseudofermion phase shifts in units of 7r given below. For the phase shifts 7r $,5^ Q „(g, q') we provide the values 
for U/t — and U/t — ► 0, when different. The expressions read, 



7r$ c o,co(q, q) = -sgn^sinfc° (g)-sinfc° (gO)^^ + ^|g|,2fc J rW[ s g I1 (9)]( 



-2V2 



1 7T . 



(C5) 



7r$ c0 ,5i(?, q ) = -sgn(sinfc° (g) - c s i(q') sin(g'))^7 



(9) 



(C6) 



n$si,co(q, q') = -sgn(sin(g) - sinfc° (g')J ^ ; q ^ ±k F 
sgn(q)n 

= 0; q = ±k F ; U/t = 0, 



(C7) 



; g^±*; f 
sgn(g)7r 



1 + 5 M /2(1-V2) ; q = ±fc F ; 17/* -► 



2^ 

= ; g = ±k F ; U/t = 0, 



(C8) 



and 



7r*c0,ci(?, ?') = -sgn(sinfc° (g) - c c i(g') sinfc^g') 



2tt 



(C9) 



7T$ s i, c i(g, q') = -0(k F - |g|)sgn sing -sin /^(g') -, 



(C10) 



35 



respectively. Here the sign function is such that sgn(O) = and 6(x) = 1 for x > and 6{x) = for x < and thus 
tt^s I, ci( d=fc-Fi = 0. In the above equations, fceo(<7) = um [//t->o where the U/t — ► value of k°(q) is given in 

Eq. ilTTlj . 



and 



(<?) = § + s gn(?) fc F ; < | Q | < [tt - 2fc F ] 
= 0; 9 = 0, 



(Cll) 



C c ( ? )=2fe(2*; F -|g|) + v^5| g | )2fcF +2fl(7r-|«|)e(|g|-2fc F )+«| 



(C12) 



c a x(q) = 1, |q|<£; F ; c s i(g) = oo, q = ±k F 

c cl {q) = 1, |q| < [7T - 2/c F ] ; c c i(g) = oo, <? = ±[7r - 2fc F ] 



(C13) 
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